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SUMMARY 

Analytical  and  numerical  studies  are  made  of  aerodynamic  forces  for  wing 
and  control-surface  motion  with  general  time-dependence  in  linearized  subsonic 
flow.  Alternative  formulations  are  discussed,  with  particular  attention  to  one 
in  the  time  domain  where  quasi-steady  displacement  and  rate  terms  are  combined 
with  a  residual  history  term.  An  accurate  calculation  procedure  is  devised,  and 
results  are  illustrated  for  a  high-aspect-ratio  wing  at  Mach  number  0.8  with 
trai ling-edge,  leading-edge  and  all-moving-tip  controls.  Emphasis  is  placed  on 
asymptotic  behaviour  at  small  and  large  times. 

Individual  control  characteristics  are  compared  over  a  wide  range  of  con¬ 
trol  rate.  The  usefulness  of  the  quasi-steady  approximation  is  established  for 
hinge  moments  and  is  analysed  for  lift,  where  the  rate  and  history  terms  become 
important  together.  The  rapid  lift  response  to  the  leading-edge  control  and  the 
sluggish  lift  response  to  the  trailing-edge  control  are  explained.  These  forces 
in  the  time  domain  are  confirmed  by  Fourier  transform  calculations  in  the  fre¬ 
quency  domain,  which  show  the  extent  to  which  the  range  of  frequency  can  be 
truncated. 

The  control-surface  motion  to  produce  a  known  time-dependent  force  is 
determined.  It  is  remarkable  how  rapidly  the  controls  can  neutralize  the  growth 
of  lift  as  the  wing  enters  a  step  gust. 
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1  INTRODUCTION 

The  concept  of  active  control  technology  has  become  associated  with  many  aspects  of 
unsteadiness  in  flight,  which  together  cover  a  wide  range  of  frequency.  The  use  of  an 
autopilot  to  control  height  and  speed  defines  the  low-frequency  end  of  the  spectrum. 

Next  comes  the  control  of  short-period  pitching  oscillations  and  the  application  to  the 
relaxed  static  stability  of  aircraft,  for  which  it  is  sufficient  to  adopt  low-frequency 
aerodynamics.  Moderate  frequencies  arise  in  the  damping  of  structural  modes  and  are 
prominent  in  aircraft  response  to  gusts,  where  active  controls  can  enhance  ride  control  or 
alleviate  structural  loads.  Higher  frequencies  may  arise  in  the  efficiency  and  safety 
aspects  of  manoeuvre  control,  and  still  higher  frequencies  are  invoked  in  problems  of 
flutter  suppression.  At  the  top  of  the  range  comes  the  fatigue  problem  of  buffet  loads. 
These  applications  make  demands  in  the  field  of  unsteady  aerodynamic  theory. 

Active  controls  are  potential ’y  of  many  types.  The  present  investigation  concerns 
wing  aerodynamics  in  attached  subsonic  flow,  and  the  treatment  covers  the  conventional 
aileron  or  trailing-edge  control  and  also  leading-edge  and  all-moving-tip  controls. 

Other  possible  types  involve  greater  theoretical  difficulty.  Discontinuities  in  planform 
arise  with  extended  flaps:  controls  such  as  spoilers  incur  problems  of  flow  separation: 
canards  or  elevators  involve  interfering  surfaces.  Whatever  control  devices  are  con¬ 
sidered,  the  technology  calls  for  their  rapid  deployment.  It  is  necessary  to  be  able  to 
calculate  the  time-dependent  aerodynamic  forces  for  arbitrary  control  motion  and,  more¬ 
over,  to  solve  the  inverse  problem  of  how  to  deploy  a  control  to  produce  a  given  unsteady 
force. 

When  linearized  subsonic  lifting-surface  theory  is  adequate,  it  is  convenient  to 

apply  the  well-established  principles  of  superposing  oscillatory  motions  to  give  step 

changes,  as  has  been  described  in  the  review  of  indicial  aerodynamics  by  Lomax1.  The 

practical  application  to  the  response  of  a  flexible  aircraft  to  gusts  is  illustrated  by 
2 

Mitchell  using  a  Fourier  transform  technique.  Subsequent  to  Refs  1  and  2  the  accuracy 

and  scope  of  numerical  methods  for  aerodynamic  lead  calculation  have  improved.  Indeed 
3 

(Corner's  recent  survey  of  theoretical  aerodynamic  methods  for  active  control  devices 
contains  over  a  hundred  references  later  than  1965  and  discusses  not  only  subsonic  flow 
but  transonic,  supersonic  and  separated  flows  as  well.  The  role  of  unsteady  aerodynamics 

,  4 

in  problems  of  aircraft  dynamics  is  discussed  by  Hancock  .  When  the  mathematical 
modelling  of  the  aerodynamics  can  satisfactorily  be  expressed  in  terms  of  the  instanta¬ 
neous  motion,  the  solution  of  the  equations  of  aircraft  dynamics  is  numerically  straight¬ 
forward.  The  synthesis  of  active  control  systems  is  the  current  problem  calling  for 
most  accuracy  and  generality  in  time-dependent  aerodynamics. 

The  present  investigation  gives  detailed  consideration  to  wing  and  control  motion 
with  general  time  dependence  in  linearized  subsonic  flow.  We  shall  look  at  alternative 
formulations  of  the  aerodynamic  forces  and  the  practical  merits  of  the  various  functions 
introduced  to  account  for  the  history  of  the  motion.  We  shall  seek  to  develop  numerical 
procedures  and  to  test  them  for  the  gross  wing  planform  and  aileron  in  Fig  1  typical  of  a 
modern  medium-sized  high-subsonic  airliner,  and  also  with  leading-edge  and  all-moving-tip 
controls.  The  characteristics  of  the  various  controls  will  be  compared  over  a  range  of 
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control  rate  for  the  Mach  number  0.8,  and  the  usefulness  of  the  low-frequency  approxima¬ 
tion  will  be  analysed.  For  a  particular  form  of  transient  control  motion,  we  shall  also 
express  the  aerodynamic  forces  as  integrals  with  respect  to  frequency,  so  that  with  the 
aid  of  a  Fourier  transform  technique  the  consistency  of  the  r.umerical  results  and  the 
permissible  truncation  of  the  infinite  range  of  frequency  can  be  examined.  As  an  illustra¬ 
tive  solution  of  the  inverse  problem,  we  shall  demonstrate  how  the  three  types  of 

control  can  be  deployed  to  neutralize  the  growth  of  lift  due  to  entry  into  a  step  gust. 

2  BASIC  TIME-DEPENDENT  WING  FORCES 

For  many  years  the  primary  requirement  for  high-frequpncy  unsteady  aerodynamics  has 
arisen  from  the  need  for  flutter  calculations,  which  continue  to  rely  heavily  on  small- 
perturbation  theories  for  wings  in  simple  harmonic  motion.  Frequency  parameter  is  there¬ 
fore  a  basic  independent  variable.  In  the  linearized  treatment  of  active  controls  one 
approach  is  to  use  the  frequency  spectrum  to  build  a  step  change,  and  hence  to  derive  a 
time  integral  to  represent  any  specified  motion.  The  corresponding  aerodynamic  forces  are 
then  obtainable  by  linear  superposition.  This  is  analogous  to  the  representation  of  an 
arbitrary  gust  profile  and  its  aerodynamics  in  terms  of  those  for  sinusoidal  gusts  via 
an  intermediate  step  gust. 

A  somewhat  different  approach  is  to  define  an  arbitrary  smooth  motion  of  a  wing  or 
control  surface  with  finite  acceleration  starting  at  time  t  *  0.  This  expresses  the 
reality  that  any  such  motion  can  be  created  by  the  application  of  a  finite  force,  which 
may  be  of  finite  duration  tj  .  Double  integral  expressions  for  resulting  aerodynamic 
forces  are  obtained  in  Appendix  A,  and  the  relative  merits  of  working  in  the  time  domain 
or  the  frequency  domain  are  discussed  in  section  2.3.  This  leads  to  the  concept  of 
hereditary  functions,  for  which  various  alternatives  are  considered  in  section  2.4. 

2. 1  Oscillatory  forces 

The  present  calculations  are  made  with  the  subsonic  lifting-surface  method  of 
Davies'3.  The  oscillatory  load  distribution  is  taken  to  be  continuous  and  finite  over 
the  planform  apart  from  the  inverse-square-root  singularity  at  the  leading  edge  and  to  be 
independent  of  the  wing  thickness  and  the  mean  flow.  For  an  oscillatory  upward  wing  dis¬ 
placement  given  by 


-  (,,,>.*«] 


the  surface  boundary  condition  imposes  a  complex  upvash  angle 


3z.  iuz. 


and  the  theory  yields  the  integral  equation 


V;(x,y) 
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Here  the  integral  is  taken  over  the  planform  where  the  typical  length  is  chosen  to  be  the 
geometric  mean  chord  c  ,  the  unknown  oscillatory  loading  over  the  planform  is  defined  by 

Lift  per  unit  area  -  fl|pU^2.j(x,y)e*utj  ,  (2-4) 

and  corresponding  to  the  circular  frequency  u  the  frequency  parameter  is 


uc 

U  * 


(2-5) 


The  main  stream  has  density  p  ,  velocity  U  parallel  to  the  x-axis  and  Mach  number  M  . 
The  kernel*  function  K  in  equation  (2-3)  is  given  by'equation  (17)  of  Ref  5  as  a 
function  of  v,  M  and  the  two  space  variables.  For  a  set  of  displacement  modes,  desig¬ 
nated  by  equation  (2-1)  with  subscript  i  ,  there  is  a  complex  matrix  of  generalized 
aerodynamic  forces  given  by  the  integrals 


Qij  "  p(/zi(VVyVVdx0dy0  (2“6) 

over  the  planform  from  tip  to  tip. 

The  displacement  or  force  modes  to  be  used  in  the  present  calculations  are  defined 
in  Table  1.  Modes  1  to  5  correspond  to  heaving  motion  and  the  four  rotational  motions  of 
the  wing  and  control  surfaces  indicated  in  Fig  1 ,  and.  equation  (2-2)  is  seen  to  apply  to 

w./U  in  Table  1.  The  theory  of  Ref  5  obtains  rigorous  approximations  to  Q..  with  the 
3  J 

replacement  of  z^.  and  w^  by  equivalent  smooth  complex  functions,  which  avoid  discon¬ 
tinuities  in  the  equivalent  at  the  hinge  line  in  modes  3  and  4.  A  modal  data 
program  to  evaluate  these  equivalent  functions  at  the  required  locations  on  the  planform 
is  described  in  section  5  of  Ref  6.  Very  simple  changes  to  that  program  make  possible  its 
extension  from  an  arbitrary  trailing-edge  control  to  one  selected  from  the  three  types  in 
Fig  1.  There  is  one  slight  complication,  however,  in  the  case  of  the  all-moving-tip 
control;  the  hinge  axis  must  not  intersect  the  leading  or  trailing  edge,  and  to  acquire 
results  for  the  forward  axis,  chosen  to  give  aerodynamic  balance,  it  becomes  necessary 
to  take  two  axes  further  aft  and  to  extrapolate  linearly  in  x^j  .  The  remaining  mode  6 
in  Table  1  is  different  from  the  other  five  in  that  w^  represents  a  sinusoidal  gust  and 
the  corresponding  quantity  Zj  that  would  be  obtained  from  the  differential  equation 
(2-2)  is  complex  and  useless  as  a  force  mode.  The  modal  data  program  is  therefore 
replaced  by  a  direct  evaluation  of  the  already  smooth  function  wj  at  the  required  loca¬ 
tions,  and  the  force  mode  is  replaced  by  the  unrelated  function  |y|/>  which  yields  a 
generalized  force  proportional  to  the  root  bending  moment. 

The  Fortran  program  for  the  numerical  solution  of  equation  (2-3)  and  the  evaluation 
of  the  generalized  force  coefficients  Q. .  in  equation  (2-6)  is  described  in  section  3  of 


*  At  a  recent  colloquium  in  Gottingen  to  celebrate  the  eightieth  birthday  of 
Professor  H.G.  Kussner,  it  was  proposed  internationally  that  this  function  should 
henceforth  be  known  as  the  Kussner  kernel  in  recognition  of  his  fundamental  work  on  the 
subject  published  in  1940. 
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Ref  6.  One  minor  change  is  made  to  both  this  and  the  modal  data  program,  so  that  the  sub 
routine  CHORD  can  handle  any  polygonal  planform  with  provision  for  arbitrary  extents  and 
alternative  shapes  of  rounding  at  the  corners;  numerical  input  data  are  read  on  the  first 
occasion  that  the  subroutine  is  entered.  The  actual  planform  is  defined  in  Fig  1,  but 
after  smoothing  around  the  central  and  trailing-edge  crank  sections  with  the  aid  of  the 
function 

g(A)  -  (1  “  A)4(5  +  4A  +  X2)  ,  0  <  A  <  1,  (2-7) 


the  modified  shape  is  as  follows.  For  the  leading  edge 

X£  v 

*  1.6934  £  +  0.1660  g(A)  , 

c 

*1  V 

—  -  1.6934  £  , 

s 
c 

with  A  *  (y/s)/0.0980.  For  the  trailing  edge 


0  <£<  0.0980 
s 


0.0980  <£  <  1 

S 


>  (2-8) 


=  1.6896  , 

c 

x, 

—  «  1.6896  +  0.0362  g(A)  , 

c 

—  -  1.4331  +  0.7380  £  +  0.0362  g(A)  , 

c 

—  -  1.4331  +  0.7380  £  , 

c 


0  <£<  0.2496 
s 

0.2496  <£  <  0.3476 
s 

>  (2-9) 


0.3476  <£  <  0.4456 
s 

0.4456  <£  <  1 
s 


with  A  -  |  (y/s)  -  0.3476|/o.0980  . 

l'n  addition  to  the  planform  data,  including  the  typical  length  c  and  the  semi¬ 
span  s  -  3.9960c  ,  it  is  necessary  to  specify  certain  parameters  to  define  the  calcula¬ 
tion.  These  control  the  output  and  ensure  that  the  number  of  modes  is  consistent  with 
the  modal  data  and  are  treated  as  synmetric  in  y  in  the  present  calculations  or  as 
antisynmetric,  if  required.  There  are  parameters  to  fix  the  numbers  of  spanwise  loading 
points  and  upwash  points,  both  of  which  are  set  to  m  *  23  in  this  instance;  similarly 
the  numbers  of  chordwise  loading  and  upwash  points  are  set  to  N  >  8.  To  enhance  the 
accuracy  of  spanwise  integration  of  equation  (2-3),  integers  q.  (i  »  1  to  N)  and  their 
lowest  common  multiple  are  increased  from  unity,  and  their  values  in  the  present  calcula¬ 
tions  are 

q.  *  (16,  8,  6,  4,  4,  6,  8,  16)  with  an  LCM  of  48  .  (2-10] 

The  Mach  number  is  set  to  0.8  and  the  frequency  parameter  takes  in  turn  the  24  values 

i 

\ 


v  ” 


0,  0.05,  0.15,  0.4,  0.6,  0.75,  1.0  and  1.2  (0.3)  6.0 


(2-11 
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Selected  force  coefficients  have  been  tabulated  against  v  .  Tables  2  to  4  give 
the  coefficients  “Qjj»  “Q2j  an<^  "Qjj  th®  trail ing-edge,  leading-edge  and  all- 

moving-tip  controls  respectively.  Each  of  these  complex  coefficients  is  split  into  its 
in-phase  and  in-quadracure  compo  »nts  written  as 


Q. . (v)  *  Q!.(v)  +  ivQV.(v) 

ij  ij  ij 


(2-12) 


The  signs  in  Tables  1  to  4  may  call  for  clarification.  The  rotational  modes  2  to  5  are 
nose-down,  as  are  the  pitching  and  hinge  moment  coefficients  with  i  *  2  and  i  •«  j  in 
equation  (2-6).  The  positive  values  of  -Qjj  represent  negative  in-phase  lift  corres¬ 
ponding  to  each  nose-down  control .  The  corresponding  in-quadrature  coefficients  -Q" j 
involve  rapid  changes  when  v  is  small,  and  for  each  control  it  changes  sign  in  the 
frequency  range  of  equation  (2-11).  The  positive  values  of  -Q^j  represent  in-phase  nose- 
up  pitching  moments  about  the  axis  x  ■  xa  in  Fig  1  corresponding  to  nose-down  control 
angles.  The  hinge  moment  coefficients  include  equal  contributions  from  the  pair  of  control 

surfaces.  They  appear  as  small  quantities  because  the  non-dimensionalizing  volume  in 
3 

equation  (2-6)  is  c  instead  of  the  more  typical  product  of  the  area  and  mean  chord  of 

3 

the  control  surfaces,  which  is  0.04c  in  modes  3  and  4.  Both  parts  of  — Qj j  change 
fairly  slowly  with  v  .  Positive  - Q indicates  a  restoring  moment  on  the  trailing-edge 
control,  unlike  -Q^  for  the  leading-edge  control.  The  all-moving-tip  is  almost 
neutral  with  static  balance  and  small  values  of  of  varying  sign.  Positive  values 

of  — Q V j  show  that  free  motion  in  any  one  of  the  three  control-surface  modes  will  decay. 


An  Argand  diagram  of  the  complex  quantity  Qj^W/Qj^O)  is  drawn  in  g  2(a) 
over  the  frequency  range  0  <  v  <  6  ,  and  it  shows  an  undulatory  behaviour  the  range 
1  <  v  <  3  and  an  asymptote  at  (^(“O/QjsW  «  0.325  ,  both  of  which  will  be  discussed 
later.  It  is  sufficient  here  to  contrast  this  behaviour  with  that  of  the  smooth  Argand 
diagram  of  Q]g('’)/Qjg(°)  corresponding  to  sinusoidal  gusts  of  varying  wavelength,  which 
spirals  to  the  origin  as  wavelength  decreases  or  v  -*■  •  . 


2.2  Indicial  aerodynamics 

The  term  'indicial  aerodynamics'  has  been  in  the  technical  vocabulary  for  decades, 
but  its  origin  is  obscure  to  many  who  use  it.  It  seems  unlikely  that  the  term  is  derived 
from  'index'  in  the  mathematical  sense,  although  the  exponent  in  equation  (2-1)  provides 
a  stepping  stone  from  single  harmonic  motion  to  one  of  growing  or  decaying  amplitude  with 
complex  frequency  u  .  It  is  probably  derived  from  'indicium'  (*  mark,  sign)  and  refers 
to  the  signature  of  a  discrete  event  in  time  such  as  a  step  change.  The  fact  remains 
that  the  term  is  applied  to  the  classical  work  of  Wagner  (1925)  and  Kussner  (1936)  who 
solved  analytically  the  growth  of  lift  in  two-dimensional  incompressible  flow  arising 
respectively  from  a  sudden  change  of  angle  of  attack  and  from  entry  into  a  uniform 
vertical  gust. 

The  developments  in  this  field  up  to  the  late  1950s  have  been  reviewed  by  Lomax'. 
An  additional  paper  of  this  period  by  Drischler^  also  gives  useful  analytical  and 
numerical  data.  Without  restriction  to  two-dimensional  or  incompressible  flow  there  are 
basic  linear  relationships  between  indicial  and  oscillatory  aerodynamics.  These  provide 


8 


alternative  expressions  for  the  lift  due  to  a  step  motion  of  a  wing  or  control  surface, 
and,  following  equations  (41)  of  Ref  I,  we  can  take  for  o  >  0 


Here  the  oscillatory  coefficients  Q.j(v)  with  i  =  1  and  j  arbitrary  are  split  into 
their  real  and  imaginary  parts  from  equation  (2-12),  and 


(2-15) 


is  the  travel  in  mean  chords  following  the  step  motion.  The  mathematical  derivation  of 
aerodynamic  forces  for  general  time-dependent  motion  and  full  details  of  their  evaluation 
are  covered  in  sections  2.3  to  3.3.  In  the  limit  as  a  -*■  »  equations  (2-13)  and  (2-14) 


pU2c2Q| j (0) 


(2-16) 


Hence  the  ratio  L(o)/L(»)  for  the  trailing-edge  control  (j  *  3)  can  be  evaluated  from 
the  lift  data  in  Table  2  or  Fig  2a,  and  the  results  are  plotted  as  the  broken  curve  in 
Fig  2b.  The  full  curve  for  entry  into  a  step  gust  is  obtained  similarly  from  equations 
(2-13)  to  (2-16)  from  the  lift  data  for  sinusoidal  gusts  (j  ■  6)  in  Fig  2a. 

Some  contrasts  and  comparisons  between  the  pairs  of  curves  for  the  oscillatory  and 
step  motions  in  Fig  2  for  M  *  0.8  will  set  the  scene.  As  will  be  explained  in  the 
following  sections,  there  is  an  impulsive  lift  at  o  ■  0  proportional  to  Q’jj(°°)  ;  this 
is  finite  and  non-zero  for  the  control  surface,  but  is  zero  for  the  gust  case  where  the 
Argand  diagram  indicates  a  curve  spiralling  to  the  origin  at  v  -*>  «®  .  The  limiting  value 
L(0)  as  o  tends  to  zero  through  positive  values  is  proportional  to  Q j j (CT>)  which, 
also,  is  finite  and  non-zero  for  the  control  and  zero  for  the  gust.  The  larger  variations 
in  curvature  of  L(o)  for  the  trailing-edge  control  in  Fig  2b  are  thought  to  follow 
from  the  distinctive  undulations  of  vQ"^ (v)  in  Fig  2a,  about  which  more  will  be  said 
later.  Similarly  it  is  worth  noting  at  this  stage  that  as  o  -*■  «®  the  two  curves  in 
Fig  2b  have  the  same  asymptotic  behaviour 


L(o)  .  1.009 

L(»)  "  2 

0 


which,  as  will  be  seen,  relates  to  a  common  feature  of  the  Argand  diagrams  near  v  *  0  . 


9 


2.3  Non-sinusoidal  motion 

In  indicial  aerodynamics  the  step  change  is  one  of  Che  usual  basic  elements.  It 
is  unrealistic,  but  serves  as  a  practical  mathematical  device  for  the  linear  construction 
of  arbitrary  time-dependent  wing  motion  or  a  gust  of  arbitrary  profile.  In  the  present 
treatment  of  non-sinusoidal  motion  the  wing  and  control  surfaces  are  regarded  as  being 
activated  by  finite  forces  and  involve  finite  acceleration.  In  the  context  of  active 
controls  it  is  therefore  postulated  Chat  the  upward  wing  displacement  in  the  mode  Zj(x,y) 
has  a  smooth  time  dependence  so  that 


z  -  qj(T)Zj(x,y) 


(2-18) 


where  x  -  Ut/c  and  (t) 
five  is  bounded. 


and  its  first  derivative  are  continuous  and  its  second  deriva- 


In  Appendix  A,  alternative  expressions  are  derived  for  the  generalized  force  in 

mode  i  due  to  the  wing  motion  in  equation  (2-18)  with  q.(x)  -  0  for  x  <  0  .  With 

J 

the  addition  of  subscripts  to  equations  (A-12)  and  (A-13),  we  have 


Q.(x) 


and 


Q^t) 


dq.(x)  ,  j  7  dq.(x  ) 

Qij(o,)  '~dx~  *  Q-j(°)qj(T)  +  -  J  |  -j--0-  iQ'ij(S)  -  QV j  (“) }  cos[v(x-x0)]dvdx 

0  0  0 


dq.(x) 

QV.(-)  — J - + 

^ij  dx 


2  j  J  dW  r-  i 

¥  J  J  - L-sin[v(x-x0)]dvdx0 

0  0  0  v 


.  (2-19) 


(2-20) 


where  Q!j(w)  and  QVj(v)  are  defined  in  equations  (2-6)  and  (2-12).  These  quantities 
can  be  calculated,  as  in  Tables  2  to  4,  for  particular  values  of  v  ,  but  they  each 
require  approximate  smooth  representation  over  the  complete  frequency  spectrum.  On  the 
other  hand  9j(tq)  ®*y  be  specified  by  a  simple  formula,  whereby  the  integrations  with 
respect  to  xQ  in  equations  (2-19)  and  (2-20)  can  be  carried  out  analytically,  so  as  to 
leave  a  single  integral  for  Q^(x)  with  respect  to  v  in  the  frequency  domain.  The  two 
double  integrals  for  Q^(x)  are  equivalent  to  those  derived  from  reciprocal  Fourier 
transform  relationships  between  the  indicial  aerodynamic  force  and  the  corresponding 
complex  coefficient  ^ (v)  fcr  oscillatory  motion  (see  equations  (5)  and  (6)  of  Ref  7). 

The  classical  approach  is  to  integrate  first  with  respect  to  v  to  give  in  the 
time  domain 


Q^t) 


dq^CO  f  d9j(T0) 

«u‘->  -fe-  *  ««<«  I  -kf-  V*- V*o  • 


(2-21) 


where  we  have  alternative  expressions 


[f] 


QV.(v)  -  QV . (“) 


Q!  .(0) 
xij 


cos  vo  dv 


(2-22) 


F,(o) 


I  [ 


—  sin  vo  dv 


(2-23) 


from  equations  (2-19)  and  (2-20)  respectively.  These  are  readily  identified  with 
equations  (2-13)  and  (2-14).  The  advantage  of  this  formulation  is  that  Fj(o)  ,  which  has 
often  been  called  a  Wagner  function,  can  be  calculated  without  reference  to  the  wing 
motion  qj(f)  .  Thus  equation  (2-21)  can  be  used  whether  qj(x)  is  regarded  as  known  or 
unknown. 

2.4  Hereditary  functions 

We  now  introduce  and  appraise  some  hereditary  functions,  which  embody  the  history 
of  the  wing  motion  in  a  general  expression  for  the  resulting  instantaneous  aerodynamic 
force.  The  function  Fj(o)  in  equation  (2-22)  or  (2-23)  is  one  of  them,  and  others  arise 
when  equation  (2-21)  is  integrated  by  parts.  Thus, 


dq.(x) 

VT>  ■  Qijw -dr-* 


T 

Qij(0)h(T0)Fi(T-V]o  -  Qij(0)  /  w 


dFi(x-y  , 

dT0  0 


dq.(x) 

-Jr- 


T 

+  Q'.(0)F,(0)q.(T)  +  Q!.(0)  J  qj(T0)F0(T-T0)dT0 


(2-24) 


where 


F0(o) 


dF , (o) 


2  f  2  f  QiiCv)  " 

■  I  j  '  qT.(O)  3 - *  sin  va  dv  »  |  j  --J  •q!;(0)  J - cos  vo  dv  (2-25) 


Q!.(v)  -  Q!.(») 


2  f  _  - 

im  —  /  — * - sin  vo  dv 

o  '5Qlj(0) 

[\  ?Q!.(v)-Q!.H  2Q ! .  («)  f  .  -  1 

.  2  I  xn'  '  I  sm  vo 

ins  —  — *■ - J -  sm  vo  dv  +  ■  - dv 


sin  vo  dv 


0  %<°> 


Q!,(~) 


P5H 


S 


3 


(2-26) 
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By  equations  (2-2£)  and  (2-26) 


Q£Ct) 


T 

qv . (“)  qj--~  +  +  f qjW’V-'o 


(2-27) 


where  FQ(o)  is  given  by  either  of  the  integrals  in  equations  (2-25). 

Another  hereditary  function  is  found  when  equation  (2-21)  is  rewritten  in  the  form 


Qj(t) 


dq.(r) 


ii<»>  *  Qii<c) I 

J  L  o  J  o  0  0 


-OdT. 


(2-28) 


where 


F2(o) 


f  Qii(W) 

J  Fl(T0)dT0  +  Q!  (0)  * 
0  1J 


(2-29) 


The  constant  in  equation  (2-29)  is  chosen  to  cancel  the  first  term  on  the  right  hand  side 
of  equation  (2-28).  Remembering  that  dqj(T0)/dT0  -  0  at  tq  -  0  ,  we  obtain 


r  <*Vv 

Q.(t)  -  Q!.(0)j  3  —  F2(T-T0)dT0 

1  l  dT„ 


(2-30) 


where,  from  equations  (2-22)  and  (2-29) 


F2(o) 


■  o  + 


00 

1/ 


QV.(v)  -  Q’.y*) 


sm  vo 


dv  + 


QV  j  (*) 

Q!.(0) 


o  + 


if 


QV.(v) 


(0)  q 


sm  vo  d- 


(2-31) 


Thus  we  have  three  hereditary  functions  Fn(o)  with  n  -  0,  1,  2,  such  that 

,2„ 


dF.(o)  d  F2(o) 
F„«  -  -a; - ^r- 


(2-32) 


The  first  of  these  from  equations  (2-25)  is  used  in  equation  (2-27)  as  a  factor  of  q^V 
in  the  integrand.  As  q.(TQ)  is  proportional  to  the  displacement  or  position  of  the 
wing  surface,  we  call  FQ(o)  the  ’hereditary  position  factor’.  Similarly,  from  equation 
(2-22)  or  (2-23),  the  ’hereditary  velocity  factor’  Fj(o)  multiplies  dqj(TQ)/dT0  in 
equation  (2-21),  and  the  ’hereditary  acceleration  factor’  F2(o)  in  equation  (2~31) 
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multiplies  d  q^T^/dig  in  equation  (2-30).  Although  the  three  formulations  are 
analytically  equivalent,  each  Fn(o)  has  its  merits  and  defects  in  numerical,  graphical 
and  practical  terms. 

The  outstanding  merit  of  the  position  factor  Fq(o)  is  that  it  applies  even  when 
q^(r)  is  discontinuous.  Equation  (2-27)  shows  clearly  the  characteristics  described  in 
the  last  paragraph  of  section  2.2,  the  impulsive  force  proportional  to  QVj(")  due  to 
the  step  change  and  the  subsequent  value  proportional  to  (®)  before  the  hereditary 
function  begins  to  contribute.  Because  Fq(o)  tends  to  zero  for  large  o  ,  the  integral 
in  equation  (2-27)  is  particularly  suitable  for  residual  effects.  However,  it  has  the 
defect  that  Fq(o)  is  very  difficult  to  compute  for  small  o  .  Moreover,  in  the  present 
application  to  active  controls  with  smoothly  varying  qj(x)  it  may  seem  unnatural  that 
the  limiting  values  Q!^(»)  and  QVj(“)  for  infinite  frequency  should  play  such  major 
roles.  However,  the  extraneous  terms  in  equation  (2-27)  can  be  transferred  inco  the 
integrand  by  adding  a  linear  coobination  of  the  delta  function  and  its  first  derivative 
to  Fq(o)  and  so  creating  a  unified  hereditary  function  (see  equation  (2-9)  of  Ref  8). 

The  velocity  function  F ^ (o)  in  equation  (2-22)  or  (2-23)  seems  to  contain  a 
satisfactory  balance  between  what  is  required  at  small  and  large  times.  Apart  from  some 
uncertainty  in  the  calculation  of  its  limiting  value  in  equation  (2-26),  it  can  be  fairly 
well  represented  graphically.  But  it  still  requires  the  addition  of  a  delta  function  at 
o  -  0  to  account  for  the  first  term  on  the  right  hand  side  of  equation  (2-21),  It  will 
be  seen  to  have  some  of  the  characteristics  of  the  original  Wagner  function,  with  the 
slight  disadvantage  that  it  remains  finite  as  o  . 

One  obvious  merit  of  the  acceleration  factor  F2(o)  is  that  it  defines  Q^(t)  in 
equation  (2-30)  without  extraneous  terms.  Moreover,  of  these  hereditary  functions  it  is 
the  easiest  to  compute  for  small  o  .  Unfortunately,  it  diverges  for  large  o  because 
of  the  first  term  on  the  right  hand  side  of  equation  (2-31).  This  difficulty  can  be  over¬ 
come  by  further  manipulation. 

If  equation  (2-31)  is  rewritten  as 
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where* 


H(o) 


00 

1/ 


Q’.'.(v)  -  Q'.'.(O)  .  - 

13  13  sm  vo 

Q  •  j  (0)  v 


(2-35) 


2  2 

This  result  presupposes  that  q(-r)  and  dq(t)/dT  vanish  at  t  *  0  and  that  d  q(t)/dT 
exists  throughout  0  <  Tq  <  t  apart  from  the  admissibility  of  discrete  discontinuities  in 
the  range.  This  last  assumption  is  compatible  with  a  pilot  applying  a  stick  force  that  is 
continuous  apart  from  isolated  sudden  changes.  The  first  two  terms  of  equation  (2-34)  are 
equivalent  to  the  force  on  a  wing  to  first  order  in  frequency  parameter  and  are  particu¬ 
larly  appropriate  while  a  control  surface  is  being  deployed  quite  slowly.  The  residual 
integral  involves  the  hereditary  function  H(o)  in  equation  (2-35),  which  we  shall  call 
the  'history  function'.  This  is  a  slight  misnomer  in  that,  although  Q^(0)qj(T) 
represents  the  steady  force  due  to  the  instantaneous  displacement  of  the  control, 

Q'.'. (0)dq. (x)/dr  includes  not  only  Che  quasi-steady  force  due  to  the  instantaneous 
J  J 

velocity  of  the  control  but  the  additional  effect  of  an  implied  low-frequency  oscillation. 
Nevertheless,  the  third  term  of  equation  (2-34)  will  be  seen  to  make  a  negligible  contri¬ 
bution  when  the  control  rate  is  fairly  small  and  slowly  varying.  In  other  situations 
H(o)  will  be  required  to  higher  accuracy  than  the  position  factor  F.(o), 

3  CALCULATION  OF  HEREDITARY  FUNCTIONS 

The  mathematical  principles  underlying  the  whole  of  section  2  are  well-established, 
but  there  is  a  tendency  to  take  their  numerical  implementation  for  granted.  Our  objective 
is  therefore  to  tackle  the  calculations  as  thoroughly  as  circumstances  permit  and  to  dis¬ 
cover  the  difficulties  that  are  likely  to  matter  in  practical  applications  of  active  con¬ 
trol  technology.  It  is  necessary  to  discover  when  approximations  may  be  expected  to 
succeed  and  when  they  are  likely  to  fail.  At  the  other  extreme  it  is  desirable  to 
establish  the  relevance  of  the  refinements  that  can  be  made. 

Practical  numerical  approximations  to  indicial  aerodynamics  have  been  developed  for 
two-dimensional  subsonic  compressible  flows  ’  ,  but  not  enough  is  known  about  high- 
frequency  effects  in  linearized  three-dimensional  flows.  It  is  not  merely  that  the  range 
of  frequency  in  experimental  work  is  too  restricted,  but  the  purely  theoretical  work  is 
seldom  taken  above  flutter  frequencies  with  consequent  lack  of  physical  insight  into  what 
to  expect.  An  object  lesson  was  taught  during  the  present  investigation  concerning  the 
unexpected  undulations  shown  in  the  Argand  diagram  of  lift  for  the  oscillating  trailing- 
edge  control  in  Fig  2a.  The  phenomenon  is  discussed  and  explained  in  Ref  9  in  terms  of 
even  larger  undulations  in  chordwise  and  spanwise  loadings  remote  from  the  control  and 
associated  with  the  acoustic  propagation  and  convection  of  disturbances  created  as  the 
control  oscillates.  For  the  present  calculations  the  undulatory  behaviour  reinforces  the 
need  for  closely  spaced  values  of  v  to  define  the  quantities  Q!j(v)  and  QVj(v)  . 

It  is  considered  that  the  24  values  of  v  in  equation  (2-11)  will  provide  an  adequate 
cubic  spline  representation. 


*  H(o)  should  not  be  confused  with  the  Heaviside  step  function. 
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In  regions  of  both  very  low  and  very  high  frequency  a  cubic  spline  fit  contradicts 
the  forms  of  the  asymptotic  expansions,  and  the  special  treatment  of  these  regions  is 
covered  in  sections  3.1  and  3.2  respectively.  The  procedures  for  calculating  the 
hereditary  functions  and  the  numerical  results  are  then  discussed  in  sections  3.3  and  3. A 
respectively.  Ideas  for  empirical  correction  factors  are  put  forward  in  section  3.5. 


3.1  Logarithmic  terms  at  low  frequency 

g 

Woodcock  has  made  an  independent  study  of  the  modelling  of  linearized  three- 
dimensional  subsonic  time-dependent  aerodynamics.  While  he  is  concerned  more  especially 
with  approximations  to  Q.j(v)  suitable  for  analysis  by  Fourier  transform  methods,  he 
also  considers  the  complex  quantity  as  a  double  power  series  in  iv  and  log(iv)  .  His 
analysis  confirms  the  leading  terms  of  the  series  to  be  used  in  the  present  investigation 


and 


2  _2  3 

Q! j (v)  m  Bq  +  B J v~  log  v  +  Bjv“  +  B^v 

_2  _2 

Q'.'.(v)  *  B"  +  B'.'v  +  Bl'v  log  v  +  B','v  , 

ij  0  12  3 


(3-1) 


(3-2) 


where  log  denotes  the  natural  logarithm.  Of  course  B'  and  B"  are  identified  with 

0  0  10 

the  known  quantities  Q!j(0)  and  Q'.'j(O)  .  Milne  and  the  present  author  have  shown 
that  both  B|  and  By  are  obtainable  from  steady  aerodynamic  forces.  Equations  (10), 
(11)  and  (16)  of  Ref  10  give 


BV 


2  b; 


iVQIj(°)Qi2(°) 


(3-3) 


the  aspect  ratio  factor  has  disappeared  because  the  reference  area  in  equation  (2-6)  is 
? 

c  in  place  of  planform  area  as  used  in  Ref  10,  and  the  change  of  sign  occurs  because 
the  upwash  mode  w^/U  in  Table  I  is  of  opposite  sign  to  equation  (15)  of  Ref  10.  The 
values  of  -Bj  are  included  in  Tables  2  to  A. 

The  extension  of  low-frequency  theory  in  Ref  10  supplies  the  leading  transient  term 
of  the  asymptotic  expansion  for  large  time  of  aerodynamic  forces  due  to  a  step  motion  in 
mode  j  .  From  equation  (26)  of  Ref  10  the  generalized  force  in  mode  i  becomes 


Qj/o) 


Qi.(°) 


q;.(o)q!2(o) 

8nc2 


+  0(c"3) 


(3-A) 


in  the  present  notation, 
equation  (3-A) 


With  reference  to  Fig  2b,  we  obtain  for  lift  with 


L(o) 

L(“>) 


Q,(o) 


q;,(o) 

n 


♦  0(o"3) 


in 


(3-5) 


The  leading  transient  term  is  independent  of  the  upwash  mode  j  ,  and  equation  (2-17) 
follows  after  the  substitution  Qj2(0)  M  -25.368  from  Table  2  of  Ref  9.  Furthermore, 
from  equations  (3-1)  to  (3-3)  and  with  reference  to  Fig  2a,  the  second  terms  in  both 
Q{ j(v)/Qjj (0)  and  vQ'^ (v) /Q j ^ (0)  are  seen  to  be  independent  of  j  . 
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The  importance  of  the  second  terms  in  equations  (3-1)  and  (3-2)  grows  with  increas¬ 
ing  aspect  ratio,  so  that  rapid  variations  in  Q*’^ (~) ,  Q'^ (v)  and  Qjj(v)  are  found  in 
Tables  2  to  4  for  the  present  wing  of  aspect  ratio  8.0.  Without  Bj  and  B"  the  leading 
transient  term  in  equation  (3-5)  would  disappear.  It  is  desirable  to  incorporate  equations 
(3-1)  and  (3-2)  in  the  calculation  of  hereditary  functions.  This  is  achieved  by  choosing 
a  fairly  small  value  v-v.  ,  by  applying  a  cubic  spline  fit  in  a  finite  region  <  v  < 
vu  ,  say,  and  by  determining  B^,  B^,  B!,'  and  B!j  to  give  continuity  in  Qjj  (v)  and 

Q'.'.(v)  and  their  first  derivatives  at  v  ■  v,  .  Thus  in  0  <  v  <  v.  we  take  equations 
ij  <■  *. 

(3-1)  and  (3-2)  with 


It  is  probably  more  important  to  curtail  the  low-frequency  end  of  the  cubic  spline 
for  Q!j(v)  than  for  QVj(v)  as  the  logarithm  only  enters  in  the  third  term  of  equation 
(3r2).  Although  there  is  no  reason  to  suppose  that  B"  is  any  smaller  than  Bj  ,  B,  is 
likely  to  have  less  influence  on  Fq(o)  in  equations  (2-25)  because  of  the  factor  v  in 
the  first  integrand. 


Piston  theory  at  high  frequency 


There  is  some  upper  limit  v  ,  say,  where  reliable  calculations  by  subsonic  lifting- 
surface  theories  become  impracticable.  For  the  present  results  at  high  v  it  was  found 
to  be  necessary  to  increase  the  number  of  chordwise  terms  from  N  *  6  to  8  and  sufficient 
to  economize  in  computation  time  by  decreasing  the  number  of  spanwise  terms  from  m  »  31 
to  23.  The  choice  vu  *  6  is  purely  a  matter  of  judgement. 

For  v  >  \>u  recourse  is  made  to  piston  theory.  Ashley  and  Zartarian*1  discuss 
piston  theory  primarily  for  its  use  in  supersonic  flow.  As  Mach  number  decreases,  the 
theory  becomes  restricted  to  increasingly  high  values  of  the  frequency  parameter,  but  it 
has  some  application  as  v  tends  to  infinity  even  in  subsonic  compressible  flow.  First- 
order  piston  theory  relates  the  upwash  angle  and  wing  loading  in  equations  (2-2)  and 
(2-4)  by  the  equation 


0  /3z.  i(i)z.\ 

-iU*v)  - 


(3-10) 


i 
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ao  that  by  equations  (2-6)  and  (2-12) 


and 


QI j (“) 


Q  ijO 


2  ffzi  ffi  dXOdyO 

"  MJJ  c  3xo  "I7” 
0 

MJJ  c  c 


(3-11) 


(3-12) 


where  the  integrals  are  taken  over  the  planform.  But  the  argument  is  fallacious.  It  is 
likely  that  both  the  real  and  imaginary  parts 
inverse  powers  of  iv  ,  so  that  we  may  write 


likely  that  both  the  real  and  imaginary  parts  of  £.  have  asymptotic  expansions  in 


(  »l  *2  \  bl  b2  b3  \ 

J  \°  iv  (iv)2  )  \°  iv  (iv)2  (iv)3  / 


(3-13) 


Thus  equation  (3-11)  only  provides  the  a^  part  of  -Q!j(*)  and  excludes  the  unknown 
higher-order  contribution  denoted  by  bj  ♦  On  the  other  hand  the  single  contribution 
from  bg  suffices  for  -Q'.\  («)  ,  and  equation  (3-12)  is  correct. 

On  substituting  the  values  of  Zj/c  and  z0/c  from  Table  1  ar.d  the  planform  data 
from  equations  (2-8)  and  (2-9),  we  calculate  Q 7. (®)  analytically  for  M  *»  0.8.  The 
smoothing  terms  with  g(A)  have  the  negligible  effect  of  decreasing  — Q** j (“)  from  19.980 
to  19.967  and  have  been  ignored  in  obtaining 


-  <)•;,(“) 

-  19.980 

(19.450) 

-  qy2(-> 

-  12.940 

(11.932) 

-  Q",(«) 

-  12,940 

(13.040) 

-  <&(-> 

-  12.478 

(11.651) 

from  equation  (3-12).  The  values  in  brackens  are  those  calculated  for  the  highest 
frequency  parameter  v  *  6  by  the  lifting-surface  theory  of  Ref  5  as  detailed  in 
section  2.1.  The  edge  conditions  of  finite  non-zero  loading  implied  by  equation  (3-10) 
are  incompatible  with  lifting-surface  theory  and  inappropriate  for  finite  frequency,  and 
several  considerations  follow  from  this.  It  becomes  clearer  why  larger  numbers  of  chord- 
wise  and  spanwise  terms  are  needed  as  v  increases,  and  why  approximation  in  the  upper 
frequency  range  is  unavoidable.  The  choice  of  N  *  8  and  m  ■  23  must  inhibit  the  accuracy 
of  the  results  for  v  *  6  to  some  extent,  but  there  is  no  reason  to  expect  rapid  con¬ 
vergence  to  the  piston  theory  limit.  Nevertheless,  the  numerical  comparisons  beside 
equations  (3-14)  are  close  enough  to  encourage  the  use  of  a  matching  process  at 

v  »  v  »  6  . 
u 


All  these  difficulties  are  intensified  in  both  the  chordwise  and  spanvise  directions 
with  the  part-span  leading-edge  and  trailing-edge  control  surfaces.  Violations  of  edge 
conditions  and  continuity  of  loading  become  especially  important!  and  the  absence  of  a 
logarithmic  behaviour  in  the  loading  at  the  hinge  from  piston  theory  brings  further 


incompatibility.  For  the 

control-surface  modes 

j  *  3,  4  and  5  ve  have 

-  Qjj(“) 

-  -QVj(-) 

'  Jp/I 

> 

(xo  “  xhj)dxodyo 

"  <>2j(-> 

-  -qj2(-> 

‘  ^/J 

^x0  '  xa)(x0  ~  xhj)dx0dy0 

► 

(3-15) 

‘  QjV”) 

'  y/J 

r(X()  -  ^)\dy0 

where  the  integrals  are  taken  over  the  combined  planform  area  of  the  port  and  starboard 
control  surfaces.  For  M  »  0.8  and  the  geometry  defined  in  Fig  1  we  calculate  analytically 
the  values  of  Q7j(®)  included  in  Tables  2  to  4.  The  following  comparisons  of  lift 


V 

~Qj3(v) 

-Q*,*4(v) 

-Q','5(v) 

6 

0.0177 

0.0309 

0.2442 

00 

0.0514 

-0.0514 

0.2331 

show  large  differences  for  the  trailing-edge  control  (j  *  3)  and  more  especially  the 
leading-edge  control  (j  *  4) .  The  case  of  the  all-moving  tip  (j  »  5)  without  any  dis¬ 
continuity  at  the  hinge  gives  a  percentage  change  from  v  *  6  to  ®  comparable  with  those 
found  for  the  wing  modes  in  equations  (3-14). 


On  the  basis  of  equation  (3-13)  with  real  coefficients  an  and  bn  ,  we  postulate 
the  asymptotic  expansions  with  leading  terms 


and 


Q!.(v) 


QVj(v) 


A! 


A2 


A0  +  17  + 

v  v 


AV  A” 

Si  *  3  *  ■£ 

V  V 


(3-16) 


(3-17) 


where  A^  is  unknown  and  Aq  equals  QVj(“)  from  equation  (3-12).  Although  Aq  cannot 
be  calculated  directly,  we  can  invoke  the  identity  in  equation  (A-10)  of  Appendix  A  as 
a  -*■  0  .  Since  the  second  term  becomes 


Q' .(v)  -  Q!.(») 


sin 


vo  dv  — -iJ- 


QM-> 


sin  vo 


dv  -*•  -  JQ! .  («) 
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the  identity  yields 


A0 


QijC) 


Q!.(0) 


Q".(-)]dv 


(3-18) 


Unless  this  formula  is  used,  equation  (3-16)  is  ineffective.  The  values  of  QL(<»)  in 

Tables  2  to  A  are  only  approxiaate  and  are  the  best  that  we  can  obtain  from  equation 

(3-18) *  but  the  results  confirm  the  fallacy  of  equation  (3-11).  Then  both  equations 

(3-16)  and  (3-17)  serve  to  match  a  cubic  spline  fit  at  the  upper  end  of  the  region 

v.  <  v  <  v  .  Continuity  in  Q!.(v)  and  QV.(v)  and  their  first  derivatives  at  v  «  v 
1  u  ij  ij  u 

is  satisfied  with 


(3-19) 


(3-20) 


(3-21) 


(3-22) 


Because  of  the  uncertainty  in  Aq  ,  implied  by  equation  (3-18),  it  seems  preferable 
to  calculate  the  hereditary  functions  from  expressions  involving  QVj(v)  rather  than 
Q!j(v)  .  This  favours  the  first  integral  of  equations  (2-25)  and  equations  (2-22)  and 
(2-31)  for  Fq(o)  with  n  *  0,  1  and  2  respectively.  In  particular  for  Fq(o)  it  should 
be  noted  that  despite  the  factor  v  the  first  integrand  of  equations  (2-25)  tends  to  zero 
as  v  tends  to  infinity. 


3.3  Numerical  procedures 

Our  next  consideration  is  the  evaluation  of  the  hereditary  functions  from  equations 
(2-22),  (2-25),  (2-31)  and  (2-35)  in  the  forms 


F0(o) 


QV.(v)  -  QV.(-) 


-i 
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2  r  QV,(v)  -qv.(-)  V 

Fl(o)  -  t+7j-  QTTT5T — co,vadv  f 


—  v  sin  vo  dv 


(3-23) 
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(3-23) 


with  the  range  of  integration  subdivided  into  0  <  v  <  v^,  <  v  <  \»u  and  vu  <  v  <  «•  . 

Different  integration  procedures  are  used  in  the  three  regions,  Simpson's  Rule  in 
0  <  v  <  v  ,  Filon's  Rule  in  v«  <  v  <  v  and  analytical  treatment  in  v  <  v  <  «  . 

I  i  U  J  U 

The  necessary  equations  are  set  out  in  Appendix  B. 

The  greatest  computational  effort  lies  in  the  evaluation  of  a  sufficient  number  of 

values  QVj(v)  by  the  lifting-surface  method  of  Refs  5  and  6.  As  found  in  Ref  9  and 

given  in  Table  2  for  the  oscillatory  trailing-edge  control  (j  *  3)  at  M  *  0.8,  it  is 

necessary  to  take  intervals  in  v  as  small  as  0.3  to  define  the  undulating  curves  against 

frequency  parameter  to  reasonable  accuracy,  especially  in  the  range  1  <  v  <  3  .  The 

quantity  v.  must  not  be  less  than  the  lowest  non-zero  value  of  v  for  which  QV.(v) 

*  ^  iJ 

is  calculated,  and  likewise  \>u  must  net  exceed  the  highest  value.  The  region  4  v  < 

vy  is  treated  by  means  of  a  standard  program  for  the  cubic  spline  through  all  the  data 

points,  excluding  v  «  0  ,  and  this  provides  interpolated  values  of  QVj(v)  at  each 

integration  point  and  at  v  «*  and  v  »  vu  the  values  of  its  first  derivative  as  well. 

To  carry  out  the  integration  procedure  for  0  <  v  <  ,  it  is  necessary  to  use  the 

low-frequcncy  solutions  for  Bq  »  Q'.\.(0)  ,  the  non-dimensional  lift  Qj ^  CO)  and  the 
coefficient  Q'.^O)  f°r  pitching  motion  so  as  to  evaluate  B**  from  equation  (3-3).  Then, 
using  Q'.'.(v)  and  its  first  derivative  at  v  ■  v  ,  we  calculate  B"  and  B"  from 

XJ  X  ^  w 

equations  (3-8)  and  (3-9)  to  complete  the  definition  of  QV.(v)  in  equation  (3-2). 
Typically  we  take  »  0.08  and  subdivide  the  range  into  16  equal  parts,  and  the  inte¬ 
gration  in  equation  (B-4)  is  a  matter  of  routine. 

For  the  region  <  v  <  vy  ,  typically  148  equal  intervals  between  0.08  and  6.00, 
the  cubic  spline  and  equations  (B-5)  to  (B-8)  are  used  to  evaluate  the  integrals.  The 
velocity  and  acceleration  factors  Fj(o)  and  F2(o)  are  served  by  equations  (B-6)  and 
(B-7)  respectively,  but  the  contribution  to  the  position  factor  Fg(o)  requires  an 
equation  like  (B-7)  with  the  frequency  parameters  v  ,  v  ,  (v  +  2kh)  and  (v  +  (2k  —  1 ) h) 
transferred  from  the  denominator  to  the  numerator.  This  procedure  with  Filon's  Rule 
involves  no  restriction  on  the  magnitude  of  o  . 

The  final  contributions  from  v  <  v  <  »  involve  07 .  (v)  and  its  first  derivative 

u  1J 

at  v-vu  and  piston  theory  for  Aq  *  QVj(“)  f*°®  equation  (3-12).  The  quantities  A" 
and  A2  are  then  calculated  from  equations  (3-21)  and  (3-22)  to  complete  the  definition 
of  QVj(v)  in  equation  (3-17).  Then  for  F^o'  to  F2(o)  the  respective  equations 
(B-ll)  to  (B-13)  can  be  evaluated  with  the  aid  of  the  integral  S^(vuo)  defined  in 
equation  (B-10).  Appendix  B  also  includes  highly  accurate  expressions  for  in 
equation  (B— 14)  or  (B— 1 5)  according  to  whether  \>uo  is  less  than  or  greater  than  7. 
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The  seeching  of  QVj(v)  at  v  ■  and  v  ■  is  illustrated  in  Fig  3.  A  program 

in  ICL  1900  Extended  Fortran  is  available  with  QVj (v)  and  related  quantities  as  input  and 
the  four  hereditary  function  in  equations  (3-23)  as  output. 

3.4  Discussion  of  results 

The  calculations  are  confined  to  one  planform  (Fig  i)»  one  Mach  number  M  -  0.8,  and 
the  lift,  pitching  moment  and  hinge  moments  associated  with  pitching  motion  (j  *  2)  and 
the  deployment  of  the  three  types  of  control  surface  (j  -  3,4,5).  These  results  are 
largely  illustrative,  but  they  will  enable  us  to  discuss  the  characteristic  behaviour  of 
the  various  hereditary  functions  and  some  general  numerical  aspects  of  their  evaluation. 

Figs  4a* 4b  show  the  position  factor  Fq(o)  ,  velocity  factor  Fj(a)  and  acceleration 
factor  Fjfa)  for  the  trailing-edge  and  leading-edge  controls  respectively.  Larger  magni¬ 
tudes  of  Fq(o)  and  Fj (o)  occur  for  the  leading-edge  control  in  Fig  4b  because  Qj^(0) 
is  so  such  smaller  than  Qj.j(0)  .  These  hereditary  functions  are  considered  in  relation 
to  their  asymptotic  expansions  for  small  and  large  o  .  By  equations  (2-25)  and  (3-17) 


► 


Thus  the  curves  of  acceleration  factor  can  be  drawn  with  confidence  near  o  *  0  .  The 
limiting  values  Fj(0)  are  subject  to  a  degree  of  inaccuracy  in  the  evaluation  of 
Q!j(“)  from  equation  (3-18),  but  the  velocity  factor  is  fairly  well  represented  near 
o*0  in  Figs  4a*4b.  Although  a  value  of  Ay  is  implied  in  the  matching  process  in  the 
lower  diagram  of  Fig  3,  for  example,  it  would  be  rash  to  quantify  it.  Therefore  the 
curves  of  Fq(o)  are  speculative  near  o  ■  0  ,  and  so  to  a  lesser  extent  are  the  undula¬ 
tions  in  the  region  o  <  3  .  Beyond  this  the  influence  of  the  region  vu  <  v  <  ® 
weakens  and  we  can  turn  to  the  behaviour  of  the  hereditary  functions  for  large  o  . 
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In  conjunction  with  the  series  in  equation  (3-2)  for  small  frequency,  equation  (2-31) 
may  be  rewritten  as 


F,(o) 


2  f  B0  + 
■a  J  q!  . 


B"  +  B'.'v  .  - 

0  1  sin  vo 


-4/S 


QV:(v)  -  B"  -  B'.'v  .  - 


QJTTOT 


LiiEJii  dv 


gi*  2bh 

0  +  QTTOT  +  nQ*'.  (0>a  +  0<O"2)  * 

ij  ij 


(3-28) 


With  Bjj  ■  QVj(O)  and  B"  from  equation  (3-3),  we  have  the  asymptotic  expansion  for 
large  o 


QV.(O)  Q!.(0)Q!  (0) 

F2(o)  “  g  +  offer  ~  8J^T-.-(0To-  + 0(q  5  » 


(3-29) 


and  by  equations  (2-32) 
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(3-30) 


F0(o) 


Q*h(0)Q!2(0)  _a 
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4«Q!j  (Oo-* 


(3-31) 


Thus  Fq(o)  tends  to  zero  the  most  rapidly.  Equation  (3-30)  is  consistent  with  the 
rigorous  analysis  of  Ref  10  leading  to  equation  (3-4).  Figs  4a  &  4b  show  how  the  charac¬ 
teristics  of  the  leading-edge  and  trailing-edge  controls,  that  are  so  different  for 
small  o  ,  become  remarkably  similar  for  large  o  .  This  is  particularly  the  case  for 
lift  when  Q!j(0)  -  Q!j(0)  and  equations  (3-30)  and  (3-31)  become  independent  of  the 
mode  j  ,  as  we  anticipated  in  equation  (2-17). 

The  remarks  in  the  paragraphs  following  equation  (2-32)  are  brought  into  perspective 
by  Figs  4a&4b.  While  for  the  trailing-edge  control  Fj(o)  in  Fig  4a  is  reminiscent  of 
the  original  Wagner  function,  the  corresponding  function  in  Fig  4b  for  the  leading-edge 
control  is  less  so.  The  lift  characteristics  of  Fj(o)  for  the  various  types  of  motion 
are  compared  in  Fig  5.  It  will  be  seen  that  for  the  pitching  motion  Fj(o)  rises 
rapidly  from  0.44  to  0.84  in  the  first  half  mean  chord  of  travel  and  thereafter  increases 
relatively  slowly.  Like  the  leading-edge  control,  the  all-moving-tip  control  gives  a 
minimum  in  Fj(o)  in  the  range  a  <  1  .  All  four  velocity  factors  for  lift  have  the 
common  asymptote  for  large  a  . 

The  remaining  hereditary  function  from  equations  (3-23)  is  the  history  function 
derived  at  the  end  of  section  2.4.  It  is  readily  evaluated  from  F„(o)  ,  and  it  follows 
that  from  equation  (3-27)  for  small  a 


Q! . (0)H(o)  ~  [QV . (°>)  -  QV.(0)]  +  [Q*  (»)  -  Q!.(0)]o  -  iA’.'o2  , 


(3-32) 
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and  that  from  equation  (3-29)  for  large  a 


H(6) 


q;.(0)Q!2<0> 

8*Q!j  (0)o 


+  0(o"2) 


(3-33) 


Typical  curves  of  H(o)  are  shown  in  Figs  6  and  7.  The  history  functions  for  lift  and 
their  comoon  asymptote  in  Fig  6  have  similar  falling  trends  in  the  range  o  >  0.5  . 

The  leading-edge  control  produces  the  most  distinctive  behaviour,  and  it  may  be  noted  that 
in  the  important  range  o  <  7  the  curve  for  the  all-moving  tip  lies  between  those  for  the 
leading-edge  and  trailing-edge  controls.  Fig  7  gives  H(o)  for  the  lift,  pitching 
moment  and  hinge  moment  due  to  the  trailing-edge  control.  It  is  worth  pointing  out  from 
Ref  10  that  the  wing  loading  associated  with  the  leading  term  in  equation  (3-33) 
corresponds  to  Q!^(0)  and  is  that  for  an  incremental  change  in  angle  of  attack.  This 
loading,  which  is  large  near  the  leading-edge  and  small  near  the  trailing-edge,  plays  an 
interesting  role  in  the  behaviour  of  the  history  functions  for  moderately  large  a  ,  so 
that  the  residual  quantities  at  a  ■  6  in  Fig  7  show  relatively  large  lift,  intermediate 
pitching  moment  and  very  small  hinge  moment. 

Next  we  consider  in  Tables  5  and  6  the  sensitivity  of  the  position  factor  FQ(a) 
and  the  history  function  H(o)  to  the  manner  of  calculation.  The  fullest  calculations 
of  lift  due  to  the  trailing-edge  control  appear  in  columns  (d) .  As  described  in  Ref  9, 
it  became  necessary  to  amplify  the  set  of  14  frequency  parameters 

v  -  0,  0.05,  0.15,  0.A,  0.75,  1.2  (0.6)  6.0  (3-34) 


used  in  the  preliminary  calculations  by  including  10  extra  values  to  give  the  set  of  24 
values  in  equation  (2-11).  Results  obtained  in  columns  (c)  from  the  14  values  of  Q'j^v) 
corresponding  to  equation  (3-34)  and  the  quantity  Q'^O")  from  piston  theory  are 
reasonably  close  to  those  in  columns  (d)  for  the  most  part.  Large  discrepancies  in  Fq(o) 
arise  in  Table  5  for  small  a  ,  but  in  Table  6  H(o)  shows  nothing  worse  than  differences 
of  0.004  in  opposite  senses  at  o  ■  2.6  and  5.0  ;  these  amount  to  about  1%,  where  dis¬ 
crepancies  of  similar  magnitude  in  Fq(o)  represent  more  than  15%.  Further  degradation 
of  the  calculations  is  found  when  the  values  of  Qj'^v)  from  lifting-surface  theory  are 
ignored  where  v  exceeds  1.2.  Columns  (a)  are  obtained  when  the  integrals  in  equations 
(3-23)  are  truncated  at  v  ■  1.2  ,  and  the  comparisons  with  columns  (d)  are  poor,  but 

within  0.01  for  a  >  6.0  ;  even  for  H(o)  in  Table  6  there  are  discrepancies  of  0.035  at 

o  m  0.6  and  0.8,  while  for  Fq(c)  column  (a)  of  Table  5  is  quite  hopeless.  Substantial 
improvement  on  columns  (a)  is  found  in  columns  (b)  when  the  value  Qj* 3  (*°)  from  piston 
theory  is  used  and  the  calculation  is  made  with  v  ■  1.2  ,  but  the  need  for  lifting- 

surface  values  with  v  >  1.2  is  still  apparent.  The  final  columns  (e)  are  obtained  when 

the  low-frequency  logarithmic  behaviour  in  equation  (3-2)  is  suppressed  by  falsely  setting 
the  quantities  Bj  and  B"  of  equation  (3-3)  to  zero.  The  position  factor  is  virtually 
unaffected  in  that  the  worst  discrepancy  between  columns  (d)  and  (e)  of  Table  5  is  0.0002. 
However,  in  Table  6  the  inaccuracy  in  H(a)  reaches  the  third  decimal  place  for 
0  >  10.0;  but,  when  it  is  borne  in  mind  that  the  setting  of  to  zero  would  not 
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suppress  Che  second  term  on  Che  righc  hand  side  of  equacion  (3-2),  Che  conclusion  is 
reached  chac  che  refinement  in  chac  equacion  is  of  liccle  praccical  significance. 

A  similar  conclusion  may  noC  apply  co  equacion  (3-1),  if  calculacions  of  Fj(o)  are 
made  from  equacion  (2-23),  This  alcernacive  Co  equacion  (2-22)  is  less  satisfaccory  on 
numerical  grounds  as  being  more  sensicive  boch  Co  Che  logarichmic  behaviour  ac  small  v 
and  Co  Che  uncercaincy  ac  high  v  implicic  in  equacion  (3-18).  Nevercheless,  a  few  com¬ 
parisons  from  Che  alcernacive  formulaeion  in  cerms  of  Q! j (v)  have  been  obcained  as 
accuracely  as  possible.  In  parcicular  calculacions  of  life  due  Co  che  crailing-edge 
control,  che  worse  discrepancies  from  equacion  (2-23)  for  Fj(o)  are  roughly  0.02  near 
o  *  0,1  and  0.01  near  o  »  1.0  ,  which  are  of  liccle  consequence  in  Fig  4a;  very  similar 
discrepancies  are  found  in  picching  mooenc  and  hinge  momenC. 

Of  che  four  heredicary  funecions,  Fq(o)  and  H(o)  in  equarions  (3-31)  and  (3-33) 

Cend  asymptocically  co  zero  for  large  o  .  The  slower  convergence  of  che  hisCory  funecion, 
-1-3 

proporcional  Co  o  inscead  of  o  for  Fq(o)  >  ref leers  Che  more  enduring  cransienc 

effects  of  an  acceleracion.  Thus  che  hisCory  eerm,  given  by  che  incegral  in  equacion 

(2-34),  makes  a  subscancial  concribucion  from  che  more  discanc  pasc  when  che  variable 

is  fairly  small.  Ic  is  a  mistake  to  suppose  chat  equation  (2-34)  can  be  integrated  by 

parts  Co  give  a  series  involving  Che  higher  derivatives  of  q.(r)  and  a  residual  incegral 

J  -1 

of  decreasing  significance.  The  behaviour  of  H(o)  proporcional  to  o  does  not  permit 
further  integration  without  introducing  a  divergent  heredicary  function  with  a  term  in 
log  a  ,  This  matter  is  discussed  further  in  section  4.2. 

3.5  Empirical  correction  factors 

The  preceding  calculacions  of  heredicary  funecions  are  linearized  not  only  in  the 
sense  that  the  basic  equations  of  Appendix  A  involve  the  principle  of  superposition,  but 
since  che  oscillatory  aerodynamics  cake  no  account  of  wing  thickness  or  viscosity.  If  wc 
suppose  that  the  key  equations  (.'-12)  and  (A-13)  can  be  applied  to  compressible  viscous 
flow  past  a  thick  wing,  empirical  correction  becomes  feasible.  The  history  function 
H(o)  is  the  most  suitable  hereditary  function  for  this  purpose,  because  of  its  smooth 
behaviour  and  the  extent  to  which  this  is  governed  by  the  asymptotic  forms  for  small  and 
large  a  in  equations  (3-32)  and  (3-33) .  Thus 


(3-35) 


The  last  of  these  relationships  implies  a  bold  assumption  that  equation  (3-3)  for  the 
first  logarithmic  term  at  small  frequency  remains  valid  in  viscous  flow. 
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1C  is  reasonable  to  suppose  Chat  experimental  or  estimated  values  of  the  steady-flow 
quantities  ^1^(0),  Qjj(0)  *nd  are  obtainable.  The  low-frequency  in-quadrature 

component  QVj(O)  may  be  more  difficult  to  estimate;  one  empirical  approach  is  that  used 
in  section  7  of  Ref  12  for  control  hinge  moments,  but  which  can  be  applied  generally.  In 
the  perpetual  absence  of  experimental  data  on  Qjj(“)  and  QVj(«)  »  it  is  probably  best 
to  regard  these  as  being  independent  of  wing  thickness  and  viscosity.  Let  us  suppose  then 
that  equations  (3-35)  are  used  to  obtain  both  theoretical  and  empirical  values  of  H(0), 
(dH(o)/da]  at  o  ■  0  ,  and  oH(o)  for  a  large.  Then  we  seek  an  empirical  correction 
factor  k(o)  to  apply  to  H(o)  . 

Clearly  the  ratios  of  the  empirical  and  theoretical  quantities  in  the  first  and 
third  of  equations  (3-35)  determine  *c(0)  and  k(<»)  ,  although  they  must  both  be  regarded 
as  somewhat  tentative.  Where  H(o)  is  a  monotonic  function  falling  smoothly  from  H(0) 
to  zero,  a  simple  practical  scheme  is  to  take 

x(o)  -  x(-)  +  lx(0)  -  tc(®)]  .  (3-36) 


In  some  instances,  such  as  the  leading-edge  control  in  Fig  6  and  the  hinge  moment  in 

Fig  7,  H(o)  rises  to  e  maximum  H(o  )  at  a  *  a  ,  say.  Thus,  perhaps,  equation 

is  to 

(3-36)  can  be  replaced  by 


k(o)  *  k(0)  , 

k(o)  «  k(<»)  +  [k(0)  -  tc(®)] 


0  <  0  v<  o 

m 

o  >  0 

m 


. (3-37) 


The  second  of  equations  (3-35)  is  available  to  check  or  refine  equation  (3-36)  or  (3-37), 
if  desired.  The  application  of  the  empirical  history  function  k(o)H(o)  would  be  to 
calculation*  .s  iciated  with  equation  (2-34) ,  where  the  experimental  or  estimated 
quantities  ‘  0)  and  QVj(O)  are  needed  in  their  own  right. 

There  is  a  simpler  alternative,  which  involves  an  empirical  value  of  lQlj(0)le 
only.  With  reference  to  equations  (2-25)  for  the  position  factor  Fq(o)  ,  we  consider 
the  assumption  that  a  constant  correction  factor 


lQii(0)1e  "  Qii(o° 
Q  r.  CO)  -  q ! .  (») 


(3-38) 


may  be  applied  to  both 


Q •  j (v)  -QVj(~) 


and  Q! .  (v)  -  Q! ^  (®) 


Then  the  factor 


K 


Q!  .(0)k 

ij 

Tqt^toji; 


1  + 


Q!j(”)  UQ£j(0)]e  "  Q*j(°)} 

lQijT0)JelQij(0)  "  ^ij^T 


(3-39) 


may  be  applied  to  all  values  of  Fq(o)  .  With  the  empirical  correction  in  this  form, 
equation  (2-27)  has  the  advantage  that  its  first  two  terms  are  unchanged.  The  same  con¬ 
stant  factor  k  would  apply  similarly  to  the  history  function  H(c)  in  equation  (2-35). 
To  use  equation  (2-34),  it  would  be  necessary  to  take  a  consistent  empirical  value 
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(3-40) 


iQVj (0)1 


e 


QV j (“)  +  <{QVj(0)  -  QV.(»)} 


QVj(O)  + 


The  meric  of  this  simpler  alternacive  would  depend  a  good  deal  on  how  well  equation 
(3-40)  fits  the  low-frequency  experimental  data. 

With  either  of  these  empirical  rules  it  would  become  possible  to  modify  the  cal¬ 
culations  in  the  following  sections.  This  would  permit  some  allowance  for  effects  that 
are  ignored  in  lifting-surface  theory. 


TRANSIENT  DEPLOYMENT  OF  CONTROLS 


A  computer  program  has  been  written  to  use  separately  the  position  factor  Fq(o) 


the  velocity  factor  Fj(o) 


and  the  history  function  H(o)  as  obtained  numerically  by 


the  procedure  described  in  section  3.3,  in  order  to  calculate  the  force  coefficient  Q.(x) 


from  the  respective  equations  (2-27),  (2-21)  and  (2-34).  The  deployment  of  a  control 

and  is  identically  zero  outside 


surface  q^(x)  is  given  algebraically  in  0  «  x  £  Xj 
this  range,  but  qj(x)  and  its  first  two  derivatives  are  chosen  to  be  continuous  for  all 
T  .  The  present  calculations  for  control  surfaces  correspond  to  the  transient  deflection 


angle 


VT) 


64 


m-rj  . 


0<  X<X 


=  0 


x  <  0  and  x  >  x 


(4-1) 


where  6^  is  the  maximum  control  angle  and  Xj  =  Utj/c  is  the  duration  of  the  excursion 
in  mean  chords  of  travel.  To  the  accuracy  of  plotting  the  calculated  quantities  Q^(x) 
from  the  three  equations  are  indistinguishable. 

Some  typical  results  for  the  trailing-edge  control  and  the  relative  contributions 
of  the  individual  terms  in  equations  (2-27)  and  (2-34)  are  discussed  in  section  4.1.  The 
influence  of  the  rate  of  deployment  or  the  duration  Xj  is  analysed  in  section  4.2. 
Section  4.3  contrasts  the  characteristics  of  the  leading-edge  and  trailing-edge  controls. 
All  these  calculations  are  in  the  time  domain,  but  section  4.4  considers  the  merits  of 
working  in  the  frequency  domain  and  the  additional  conclusions  that  can  be  drawn  from 
this. 


4.1  Forces  due  to  trailing-edge  control 

The  concept  of  active  controls  envisages  the  facility  to  deploy  them  rapidly. 

The  really  significant  range  of  Xj  is  probably  5  4  Xj  <  40  .  For  example,  in  equai  urn 
(4-1)  Xj  =  10  would  give  a  maximum  control  rate 


max 


>.qJ 


6qU 

V 


5  U 

0.343  ~  , 
c 


(4-2) 


which  is  about  100  deg/s  with  6^  =  5  deg,  U  «  240  m/s  and  c  =  4  m  .  With  Xj  =  10 
the  results  for  the  lift,  pitching  moment  and  hinge  moment  due  to  the  trailing-edge 
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control  at  K  ■  0.8  are  illustrated  in  Fig  8.  The  calculations  are  made  with  unit  6^  , 
and  the  coefficient  of  lift  Qj(t)  only  just  reaches  807.  of  its  peak  quasi-steady  value 
Qj-j(O)  .  Moreover,  this  occurs  at  time  t/tj  *=  t/Tj  =  0.535  and  lags  behind  the  control 
angle  which  peaks  at  t/Tj  *  0.5  .  There  is  a  correspondingly  slow  decay  of  the  residual 
lift  after  the  motion  has  ceased.  In  contrast  to  the  lift,  the  h^nge  moment  overshoots  its 
peak  quasi-steady  value  by  4%  at  t/tj  11  0.45  and  thus  leads  the  control  angle.  There  is 
a  corresponding  overshoot  where  Q^iO/Q^iO)  becomes  negative  before  the  motion  is 
completed,  and  the  residual  effect  after  t/Tj  =  1  is  negligible.  As  might  be  expected, 
the  behaviour  of  the  transient  pitching  moment  has  intermediate  characteristics. 

Fig  9  shows  the  individual  contributions  of  the  terms  in  equations  (2-27)  and  (2-34) 
to  the  lift  when  Tj  ■  5  and  20  respectively.  The  lift  now  reaches  respectively  72%  and 
89%  of  its  peak  quasi-steady  value,  here  denoted  by  Lq  ,  at  time  t/Tj  =  0.545  and  again 
lags  behind  the  motion.  When  the  history  function  is  used  in  equation  (2-34),  the  first 
term  denotes  the  quasi-steady  contribution  proportional  to  6  .  The  second  term  is 
proportional  to  d6/dt  and  turns  out  to  be  of  opposite  sign;  this  gives  a  qualitative 
explanation  of  the  lag  in  the  peak  lift.  However,  there  remains  an  important  contribution 
from  the  history  term,  indicated  by  the  vertical  arrows  in  the  upper  diagram  of  Fig  9, 
which  changes  sign  twice.  For  the  more  rapid  motion  the  lower  diagram  similarly  disects 
equation  (2-27)  with  the  infinite-frequency  coefficients  and  the  position  factor.  It 
can  be  seen  that  the  major  contributions  are  from  the  first  term  proportional  to  d6/dr 
when  0  <  t/Tj  <  0-95  and  from  the  second  term  proportional  to  6  when  0.05  <  t/Tj  < 
0.40.  But,  thereafter,  the  position  integral  provides  the  major  contribution  which  remains 
positive.  Apart  from  the  dominant  role  of  the  very  high  frequency  parameters  and  the 
irrelevance  of  the  low-frequency  coefficients  Q} 3(0)  and  Q'^ CO)  in  the  very  early 
stages,  the  analysis  in  terms  of  the  history  function  has  rather  more  to  commend  it. 

We  would  hope  to  find  circumstances  where  the  first  two  terms  of  equation  (2-34) 
suffice,  so  that  the  history  term  can  be  ignored.  Suppose,  for  example,  that  such 
quantities  as  those  in  Figs  8  and  9  are  required  to  an  accuracy  of  about  ±0.05,  and  con¬ 
sider  the  following  cases  at  the  extremes  of  the  significant  range  of  Tj 

(a)  Lift,  Tj  =  5, 

(b)  Lift,  Tj  =  40, 

(c)  Hinge  moment,  Tj  *  5. 

The  quantities  from  the  control-angle,  control-rate  and  history  terms  of  equation  (2-34) 
are  given  in  Table  7.  In  case  (a)  the  control-angle  term  neglects  contributions  as  high 
as  0.3;  moreover,  the  history  term  is  mostly  of  larger  magnitude  than  the  total  lift  and 
its  main  effect  is  to  cancel  large  superfluous  contributions  from  the  control-rate  term. 

For  such  a  rapid  motion  at  Tj  =  5  the  analysis  of  equation  (2-34)  for  lift  ceases  to 
be  helpful.  For  the  more  gradual  deployment  with  Tj  *  40  in  case  (b) ,  Table  7  is  more 
encouraging.  The  control-angle  term  neglects  contributions  between  ±0.1  in  lift,  while 
the  inclusion  of  the  control-rate  term  reduces  the  gaps  to  within  ±0.05  or  thereabouts. 

Under  such  conditions  the  individual  terms  of  equation  (2-34)  become  useful  and  only  a 
rough  approximation  to  the  history  term  is  needed.  The  rapidly-produced  hinge  moment  in 
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case  (c)  is  an  instance  where  the  control-angle  term  neglects  positive  and  negative  contri¬ 
butions  of  magnitudes  greater  than  0.3,  while  those  of  the  history  term  hardly  exceed 
0.05.  This  special  result  is  closely  linked  to  the  small  dependence  of  Q^-j  and  Q1^  on 
v  in  Table  2,  whereby  in  Fig  7  the  history  function  itself  becomes  small. 

4.2  Duration  of  deployment 

Before  discussing  in  more  detail  the  effect  of  Xj  on  the  present  calculations,  we 
consider  in  general  terms  the  influence  of  control  rate.  If  a  control  angle  is  changed 
slowly  from  0  to  6j  ,  it  is  reasonable  to  express  a  consequent  time-dependent  force  in 
the  form 

Q.(t)  -  Q!.(0)«(t)  +  QR^r+  Qh  (4“3) 

with  some  rate  factor  .  The  boundary  condition  representing  an  instantaneous  control 
rate,  equivalent  to  the  second  term  on  the  right  hand  side  of  equation  (2-2),  can  be 
treated  separately  from  the  first  term  as  if  it  were  a  steady  state.  This  would  determine 
a  value  for  Q_  different  from  the  low-frequency  quantity  Q'.'.(0>  in  equation  (2-34). 

K  1 J 

If  this  steady  control-rate  concept  were  used  for  ,  the  residual  term  would 
represent  the  history  of  the  motion  as  the  influence  of  time  dependence  in  the  surface 
distribution  of  upwash  angle.  Because  the  very  existence  of  a  control  rate  implies  a  time- 
dependent  upwash,  such  an  approach  is  illogical  and  will  be  shown  to  be  ineffective.  If 
the  deployment  from  0  to  6j  is  carried  out  very  rapidly,  we  have  already  established 
from  equation  (2-21)  that  the  velocity  factor  F,(o)  defines  the  instantaneous  force. 

Thus  for  x  >  0 

Q.(t)  -  Q^iOfijF,  (x)  .  (4-4) 

If  the  deployment  takes  time  tj  ■  tjC/U  ,  then  a  good  approximation  to  the  force  for 
t  >  Tj  is 

Q!.(°)6  J 1 

vt)  “  I  VT-VdT0 

o 

Q! .(0)6, 

*  — ^ -  (f2(t)  “F^x-Tj)]  •  (4-5) 

A  more  crude,  but  still  useful,  approximation  is 

Qi(T)  “  Q-j<°>6,F|(T- &V  (4-6) 

which  tends  to  the  limit  in  equation  (4-4)  as  Xj  decreases  to  zero. 

Fig  10  shows  the  time-dependent  lift  as  the  trailing-edge  control  angle  6  follows 

equation  (4-1)  for  different  values  of  tj  .  The  quasi-steady  limit  Xj  -*•  «®  identifies 

L/Lq  with  6/6q  .  The  time  scale  is  normalized  with  respect  to  the  excursion  time. 

The  main  characteristic  of  Fig  10  is  the  progressive  fall  in  maximum  L/Lq  from  unity 

to  0.72  as  x.  decreases  to  5.  Subsidiary  features  when  x,  $  40  are  the  lag  in  the  '< 

i  1  » 

i, 
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occurrence  of  suximum  L/Lg  from  y/Tj  *  0.50  to  about  0.545  and  the  greater  residual 
lift  after  the  control  angle  has  returned  to  zero.  As  supplementary  information  the 
following  table  gives  the  normalized  times  w*'en  c^e  residual  lift  has  fallen  to 

L/LQ  -  0.01  . 


T1 

80 

40 

20 

10 

5 

VT1 

1.00 

1.12 

1.36 

1.83 

2.68 

The  asymptotic  limit  for  small  Tj  corresponds  to  ^(t^)  "  0.01  ,  which  in  this  case  is 
found  to  be  *  8.5  . 

The  contributions  to  Fig  10  from  the  history  term  in  equation  (2-34)  are  plotted  as 
AL/Lq  in  Fig  II.  The  widely  varying  curves  against  t/tj  -  i/Tj  illustrate  different 
aspects  of  the  application  of  the  history  function.  Provided  that  AL/LQ  is  everywhere 
small  compared  with  unity,  the  isolated  history  term  usefully  provides  a  correction  to 
the  quasi-steady  forces,  and  this  is  the  case  for  Yj  20  .  It  is  interesting  to  observe 
how  the  curves  change  shape  from  Yj  *  5  to  Tj  *  80  . 

2  2 

When  Tj  is  small,  H(y-yq)  is  slowly  varying  compared  with  d  q.(Yg)/dYg  ; 
thus,  when  t  <  Yj  ,  the  history  term  for  lift  may  be  written  as 

..  }  ‘‘2,i<,0>  dqiW 

r0  •  H"-Vj  •  U<,'V  -h~  <4-?> 


where  tq  is  some  average  value.  Similarity  between  AL/LQ  for  Yj  *  5  and  dqj(x)/dY 
is  apparent,  and  with  further  decrease  in  Yj  the  curve  becomes  increasingly  anti¬ 
symmetric  about  t/T|  ■  0.5  and  its  maximum  and  minimum  approach  y/T|  «  0.276  and  0.724 
respectively,  those  of 


(4-8) 


from  equation  (4-1).  Moreover,  from  equation  (4-7)  these  stationary  values  will 
ultimately  behave  like 


H(0) 


rdv,r| 

QV.(-)  -  Q"(0> 

‘3.43* 

.  +  hi 

L  dT  J 

.  "  -  "  q! .(6) 

max 

L  Ti  J 

“  ti 

(4-9) 


for  i  *  1  and  j  *  3  ,  which  gives  values  ±1.04  when  Tj  *  5 
and  -I .03  in  Fig  1 1 . 

The  limiting  shape  for  large  Yj  has  academic  interest, 
that  is  slowly  varying  compared  with  H(y-Yq)  ,  so  that 


to  compare  with  +0.83 

2  2 
It  is  then  d  q^T^/dx^ 
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H(T'To)dTo 


where  Tq  is  some  average  value  and,  as  Xj  increases,  Tq/t j  approaches  x/Xj  .  As 
AL/L_  decreases  with  increasing  t,  ^  its  shape  in  Fig  II  tends  to  become  synmetrical 

u  2  2  * 

like  d  qj(x)/dx  and  its  zeros  approach  the  positions  x/x^  *  0.276  and  0.724  .  The 
convergence  is  slow  and  is  complicated  by  the  fact  that  the  integral  in  equation  (4-10) 
is  logarithmically  divergent  as  t  -*•  «  .  Nevertheless,  for  motions  of  finite  duration, 
equation  (2-34)  can  be  expressed  in  the  form 


dq.(t)  dq.(i) 

Q.(t)  *  Q!.(0)q.(T)  +  QV.(0)  — —  +  Q! . (0) J(0)  - L— 

l  ij  3  ij  dr  ij  ^  2 


where 


*  f 


J(0)  +  lH(o0)do0 


J(T-T0)dT0 


(4-11) 


(4-12) 


The  practical  significance  of  equation  (4-11)  is  slight,  because  its  only  relevance  is  to 
situations  where  the  history  term  is  already  becoming  negligible.  Nevertheless,  equations 
(4-10)  to  (4-12)  yield  the  academic  result 


dq.(x)  d2q.(x) 

Q£(t)  *  Q!j(0)q.(x)  +  QV.(0)  -j-—  +  Q!j(0)[j(x)  -  J(0)]  — 


(4-13) 


for  large  x(  . 


Fig  12  shows  how  the  curves  of  instantaneous  centre  of  pressure  due  to  the  trailing- 
edge  control  vary  with  Xj  .  There  are  three  limiting  values 


for  small  x  , 


for  large  x  ,  — 


quasi-steady. 


Q^(-) 

Q-W 

Q^2(0) 

q]^W 

q;,(o> 


(4-14) 


the  first  of  which  follows  from  the  first  term  on  the  right  hand  side  of  equation  (2-24) . 
The  second  limit  follows  from  the  history  term  of  equation  (2-34)  in  conjunction  with  the 
asymptotic  expansion  in  equation  (3-33)  and  is  identified  with  the  wing  aerodynamic 
centre.  Neither  of  these  limits  is  of  much  practical  significance,  because  in  both  cases 
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Che  associated  life  force  is  very  small.  Nor  is  the  quasi-steady  value  of  much 
consequence,  for  it  belies  the  fact  that  in  the  most  important  range  5  <  Tj  <  40  and  even 
for  Tj  ■  80  the  quantity  uf/cL  is  rapidly  varying.  It  is  also  of  interest  in  Fig  12 
that  the  limit  for  large  t  is  approached  from  below;  thus  for  Tj  ■  5  and  10  the 
residual  lift  is  acting  well  forward  of  the  aerodynamic  centre. 

Tables  8  and  9  have  been  prepared  to  develop  the  theme  of  Table  7,  discussed  at  the 
end  of  section  4.1,  and  to  consider  the  influence  of  Tj  on  the  usefulness  of  equation 
(2-34)  with  the  history  term  omitted 

dq.(x) 

Qi(T)  **  Q£j<°>9j<T)  +  Q'j j  (°)  — Jj—  .  (4-15) 

Table  8  gives  the  maximum  and  root  mean  square  (rms)  errors  from  equation  (4-15)  with  and 
without  its  second  terra  for  the  20  values  t/Tj  *  0.05  (0.05)  1.00  .  As  far  as  lift  is 
concerned,  the  addition  of  the  second  term  often  increases  both  errors,  but  clearly  it 
will  bring  improvement  if  tj  is  large  enough.  Table  9  lists  the  rms  errors  and  the 
values  of  Tj  when  equation  (4-15)  changes  from  being  more  inaccurate  to  being  more 
accurate  as  a  result  of  including  the  second  term;  it  also  gives  the  lower  bounds  of  Tj 
for  which  equation  (4-15)  with  either  one  or  two  terms  gives  rms  errors  below  5%  of  the 
peak  steady  value.  For  the  case  of  lift  due  to  pitching  motion  there  is  no  range  of  Tj 
for  which  equation  (4-15)  achieves  an  rms  error  below  57.  and  its  first  term  alone  fails  to 
do  so.  For  lift  with  trailing-edge  control  motion  there  is  such  a  range  36  <  Tj  <65  ; 
by  contrast  the  corresponding  range  for  hinge  rwient  is  4  <  Tj  <  24  . 

A  subsidiary  study  is  made  to  optimize  the  rms  error  from  equation  (4-3)  with 
Qh  M  0  and  variable  QR  .  The  following  results  are  obtained  for  the  three  examples  in 
Table  7: 


(a) 

Lift,  Tj  *  5, 

Qr" 

0.293  -  Q”  (1.0) 

(b) 

Lift,  Tj  *  40, 

Qr- 

1.913  -  Q y3  (0.15) 

(4-16) 

(c) 

Hinge  moment,  Tj  -  5, 

-0.00969  -  Q£j  (1.4)  . 

Far  from  relating  to  the  rejected  steady  control-rate  concept,  which  would  require  values 
Qr  *  -0.187  for  (a)  and  (b)  and  -0.00420  for  (c),  equations  (4-16)  state  that  the  minimum 
rms  errors,  occur  with  Q„  *  Q'.'.(v)  for  frequency  parameters  v  which  are  found  to  be 
fairly  close  to  2n/Tj  corresponding  to  an  oscillation  of  period  tj  *  TjC/U  .  Neverthe¬ 
less,  the  general  approximation 

VT>  *  Qij(^7)MT)  +  Qi’j(l7)  ^r1 

would  be  a  poor  substitute  for  the  full  calculation,  unless  the  frequency  parameter 
2n/Tj  is  predominant. 

4.3  Type  of  control 


t 


Although  the  load  distribution  due  to  each  of  the  control  surfaces  is  greatly 
influenced  by  the  wing  tip,  it  is  instructive  to  recall  the  steady  two-dimensional  load 


8iUbO 


distributions  associated  with  each  type  of  control.  In  linearized  incompressible  flow  the 

2 

lift  per  unit  area,  JpU  AC^  ,  for  a  control  angle  6  is  given  for  the  trailing-edge 
control  by 


for  the  leading-edge  control  by 


-  4’h)  cot  1$  -  log 


11*  "  *hl 

i(*  +  th) 


sin  -  4>.  | 


si 

cot  *♦ +  l08|inT(? 


-  vn 
^■jJ  * 


B 


(4-18) 


(4-19) 


and  for  the  all-moving-tip  control,  as  for  a  flat  plate,  by 


AC 

*■  4  cot  .  (4-20) 


In  equations  (4-18)  to  (4-20)  with  x  and  x^  measured  from  the  leading-edge  the  angles 
<j>  and  <j>^  are  defined  as  usual  by 


^  •  JO  -  cos  <|>) 

T  "  i(1  “  cos  ♦h)  . 


(4-21) 


2  -} 

If  required,  the  simple  compressibility  factor  (1  -M  )  1  enters  on  the  right  hand  sides 
of  equations  (4—1 8)  to  (4-20).  The  load  distributions  at  M  *  0  for  the  trailing-edge 
control  with  x^  *  0.75c  ,  the  leading-edge  control  with  x^  ■  0.25c  and  the  all-moving 
tip  are  drawn  in  Fig  13.  Positive  6  is  defined  as  nose-up,  so  that  the  discontinuities 
in  surface  slope  at  the  hinge  are  in  opposite  senses  for  the  leading-edge  and  trailing- 
edge  controls.  It  follows  from  equations  (4-18)  and  (4-19)  that  for  a  given  control 
chord,  not  only  are  the  logarithmic  singularities  at  the  hinge  equal  and  opposite,  but 
the  inverse-square-root  singularities  at  the  leading-edge  are  identical  in  the  two  cases. 
Fig  13  also  shows  that  the  loading  aft  of  the  hinge  of  the  leading-edge  control  remains 
negative.  The  cancellation  effect  is  such  that  the  two-dimensional  lift  due  to  the 
leading-edge  control  is  only  9.5%  of  that  produced  by  the  trailing-edge  control;  the 
ratio  9.3%  in  three  dimensions  for  the  particular  planform  at  M  »  0.8  is  remarkably  close. 
The  familiar  flat-plate  distribution  for  the  all-moving-tip  control  in  Fig  13  is  less 
relevant,  because  the  centre  of  pressure  near  the  tip  moves  well  forward  of  quarter  chord 
to  the  extent  that  the  hinge  line  in  Fig  1  is  taken  to  balance  the  control. 

The  calculations  of  lift  due  to  the  transient  deflection  anf,le  in  equation  (4-1)  are 
recorded  in  Table  10  when  Tj  *  5  and  40  for  pitching  motion  6/0q  and  for  deployment 
of  the  three  control  surfaces  at  H  «  0.8.  For  the  slower  motion  with  Tj  “40  ,  the 
trailing-edge  control  stands  out  as  the  slowest  in  producing  and  relinquishing  lift. 
Pitching  and  leading-edge  control  rotation  give  remarkably  similar  time  histories  of 
L/Lq  ,  while  that  of  the  all-moving  tip  is  intermediate  to  the  other  two  controls. 


This  pattern  is  largely  repeated  for  the  rapid  notion  Tj  ■  5  .  The  leading-edge  control 
shows  rapidity  of  aerodynamic  response  in'  that  the  lift  overshoots  the  quasi-steady  range 
0  <  L/Lq  4  1  .  This  behaviour  is  nearly  matched  by  the  case  of  pitching  notion  where 

L/Lq  *  -0.012  at  t  m  4.5  •  As  we  have  already  seen,  the  trailing-edge  control  produces 

a  sluggish  behaviour  in  lift,  which  only  reaches  121  of  the  quasi-steady  maximum  and 

retains  135!  of  this  at  the  conclusion  of  the  control  motion.  Figs  10,  14  and  15  show  the 

influence  of  excursion  time  on  lift  for  the  trailing-edge,  leading-edge  and  all-moving- 
tip  controls  respectively. 


The  results  for  Tj  «  5  can  be  analysed  further  in  relation  to  equation  (4-17), 
which  gives  for  unit  6q 


Q,  <T) 


Q-W 


Q|  • (2s/t  j) 
-  ■ 


T(6f-  [d(T)  +Q^TS7VjT^7 


6(t) 


Qi'jL(2Tr/T1)  d6 m 


(4-22 


For  convenience  we  shall  take  the  approximate  value  2it/i|  ■*  1.2  and  use  the  data  in 
Tables  2  to  4.  It  can  be  seen  from  the  following  table  that  equation  (4-22)  provides  a 
qualitative  description  of  the  differing  characteristics  of  the  three  controls. 


fl“5 

3-3 

TE  control 

j  -4 

LE  control 

j  -5 

All-moving-tip 

Peak  L/Lq 
at  t/tj  « 

0.72 

1.05 

0.81 

0.54 

0.40 

0.45 

Q'jO^/Qj^O) 

0.611 

0.940 

0.728 

Q?j(1.2)/Q«j(1.2) 

-0.225 

0.688 

0.274 

From 

Peak  L/Lq 

0.87 

2.75 

1.14 

equation  (4-22) 

at  t/t  j  “ 

0.64 

0.30 

0.34 

The  peak  values  of  L/Lq  from  Figs  10,  14  and  15  are  ordered  according  to  the  values  of 
the  factor  Q|j(1.2)/Q|j(0)  .  The  location  of  the  peak  along  the.  time  scale  is  ordered 
by  the  ratio  Q"j(1.2)/Qj j(1.2)  which  governs  the  influence  of  the  rate  term.  When  this 
ratio  is  negative  for  j  ■  3  the  peak  is  seen  to  be  delayed,  and  when  the  ratio  is 
larger  and  positive  for  j  ■  4  the  peak  is  advanced  from  0.50  to  a  larger  extent.  But 
the  quantitative  predictions  of  equations  (4-22)  are  shown  to  caricature  the  behaviour 
and  are  included  to  emphasize  the  particular  danger  of  using  oscillatory  aerodynamics 
for  a  discrete  frequency  in  response  studies  of  active  control  systems. 

Fig  16  contrasts  the  normalized  time-dependent  lift  from  leading-edge  and  trailing- 
edge  controls  when  deployed  very  briefly  with  Tj  ■  2  .  This  could  correspond  to  a  rapid 
on-off  situation  over  a  small  range  of  control  angle.  The  trailing-edge  control  only 
reaches  585!  of  the  quasi-steady  peak  lift,  while  the  leading-edge  control  gives  a  peak- 
to-peak  variation  that  is  nearly  three  times  the  quasi-steady  one.  Clearly  the  growth  of 
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lift  remote  from  the  control  surface  is  subject  to  the  major  time  delay.  Fig  *  it.dicti'.fts 
that  with  a  positive  control  rate  this  delay  gives  lift  reduction  for  the  trailinfrad£fc 
control  and  ?ift  enhancement  for  the  leading-edge  control,  and  the  effect  is  revrrsed  as 
the  controls  return  to  their  initial  positions.  Similar  reasoning  leads  us  to  expect  the 
fairly  correlation  between  time-dependent  hinge  moments  and  their  quasi-sU**dy 

values  implied  in  Tables  8  and  9.  For  the  leading-edge  control  this  correlation  is  par¬ 
ticularly  close,  because  u.  ■  ratio  Q^(v) /Q^(v)  from  Table  3  is  as  small  as  -0.06  for 
v  <*  0  end  virtually  remains  between  half  and  twice  this  value  over  the  whole  frequency 
rang.  ,  There  seems  to  be  little  point  in  isolating  the  history  terms  for  the  lift  due 
to  the  leading-edge  and  all-moving-tip  controls  because,  as  Tables  8  and  9  show,  the  two- 
tens  approximation  from  equation  (2-34)  has  virtually  no  advantage  over  the  one-term 
approximation  if  55!  rras  accuracy  is  soofHi:,  Full  calculations  are  needed  if  tj  <  30  . 
and  a  similar  conclusion  can  he  reached  for  pitching  moment  as  well  as  lift. 

4.4  Use  of  Fourier  transform 

We  now  return  to  equations  (2-19)  and  (2-20),  derived  in  Appendix  A.  Instead  of 
integrating  with  respect  to  frequency  parameter  and  using  hereditary  functions,  we  con¬ 
sider  the  evaluation  of  Q^(t)  in  the  frequency  domain,  when  q^(r)  is  defined  by 
equation  (4-1).  It  is  more  convenient,  however,  to  work  in  complex  numbers  and  to  take 


«  «  iv(r-T  ) 

r  t  dq;(0  Q. . (v)e  1 


.  r  t  dq.(T-)  Q. 

jWQyW-i//  -ir2~± 

r\  r\  v) 


dT0dv 


(4-23) 


from  equations  (A-5)  to  (A-7)  of  Appendix  A.  From  equation  (4-1)  the  analytical  integra¬ 
tion  with  respect  to  yields 

,  fd<l j(V  ^0  „ 

I(vV  “  i  dT0 

0  u 

.  64  [(_&_,  .  _*•  }  (.  .-Wt>  .  ,)  .  (.-*'■  .  ,j  . 

l(-ivTjr  (-ivTj)5]  \  /  l(-iVTj)4  (-iVTj)6J  \  t 

.  (4-24) 


If  vtj  <  1.25,  I(vtj)  is  better  evaluated  from  the  expansion 


I(VTj) 


Rti 

♦{« 


-_  ,2  4  N4  .  2  n6  1,-^8.  1 

VTP  ~  315  (VTP  7425  (vV  “  315315  (vV  41277600  (vT 


.  16  4  .3  1  ,5  1  .7  1  .9 

1  35  "  63  495  (vTP  "  32175  3439800 


-  jol 

vt,)  J 


(4-25) 
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Since  <lj(®)  ■  0  ,  equation  (4-23)  reduces  to 


«-4/ 


•  t  Q* • (v)I(VT.)  .- 

i  f  ZiJ - —  e1VT  dv 


r  iQ.  .(v)I(vx.)  . 

I  ,  — _Lt  1  ,  e1 

2irv 


1VT  dv  ,  (4-26) 


where  Q.j(-v)  is  the  complex  conjugate  of  Q.j(v)  .  Equation  (4-26)  states  that  Q.(t) 
is  the  inverse  Fourier  transform  of  the  expression  in  curly  brackets. 

13 

A  fast  Fourier  transform  technique,  based  on  the  work  of  Cooley  and  Tukey  ,  has 
been  used  to  evaluate  equation  (4-26)  with  the  aid  of  a  cubic  spline  to  represent  the 
complex  generalized  force  coefficient  Q.j(v)  and  equations  (4-24)  and  (4-25)  for 
I (vt j )  .  In  these  calculations  it  is  necessary  to  truncate  the  range  of  v  and  straight¬ 
forward  to  vary  the  upper  limit  vu  .  The  lower  the  value  of  T|  ,  the  higher  vy  must 
be.  With  Tj  ■  10  and  \>u*5.89  the  full  calculations  of  lift  due  to  the  trailing-edge 
control  from  Fig  8  are  reproduced  to  an  accuracy  of  ±0.001.  With  t  j  “  5  near  its 
practical  lower  limit,  results  are  given  in  Table  11  for  four  values  of  vy  in  the  range 
1.23  4  \>u  i  4.91  .  For  vu»4.91  the  accuracy  of  the  Fourier  transform  calculations  is 
±0.0012  and  the  rms  error  in  the  range  Q.5  <  t  <  5.0  is  0.1%.  Table  11  shows  how  this 
rms  error  increases  to  9.6%  as  \>u  is  reduced  to  1.23.  To  keep  the  errors  within  ±0.01 
of  the  maximum  quasi-steady  lift,  it  is  advisable  to  keep  \>u  to  2.4  or  above.  For  most 
practical  purposes,  therefore,  the  range  of  frequency  parameter  0  ^  v  «  2.4  should 
suffice.  The  upper  limit  is  roughly  proportional  to  1/tj  ,  so  that  for  many  applica¬ 
tions  it  would  only  be  necessary  to  consider  v  up  to  flutter  frequencies;  on  the  other 
hand,  the  extreme  examples  in  Fig  16  with  Tj  *2  require  the  whole  range  of  the  lifting- 
surface  calculations  in  Tables  2  and  3. 

Unless  the  calculations  involve  step  changes  in  control  displacement  or  rate,  the 
oscillatory  data  for  infinite  frequency  are  clearly  irrelevant  in  practical  cases*.  It 
follows  from  equations  (3-35)  that  the  behaviour  of  the  history  function  H(o)  as  o  -*■ 
is  likewise  irrelevant.  Similar  remarks  apply  to  the  position  and  velocity  functions 
Fq(o)  and  Fj(o)  with  greater  force,  and  with  particular  reference  to  equation  (3-24) 
no  serious  concern  need  be  felt  about  the  inability  to  produce  reliable  curves  of  Fq(o) 
in  Figs  4a&4b.  Some  idea  of  the  importance  of  Fj(o)  for  small  a  can  be  gleaned  from 
equation  (4-6)  by  setting  a  suitable  lower  limit  to  the  time  required  to  deflect  a  control 
surface. 

Latitude  in  the  accuracy  to  which  high-frequency  aerodynamic  data  are  required 

g 

encourages  alternative  approaches.  Thus  the  aerodynamic  modelling  discussed  by  Woodcock  , 

which  exploits  Fourier  transforms,  can  be  seen  in  a  practical  light.  The  analysis  that 

follows  from  representation  of  the  coefficients  Q..(v)  as  the  ratio  of  polynomials  in 

14 

v  ,  as  proposed  by  Vepa  for  example,  becomes  credible.  Although  such  an  approximation 
could  hardly  be  expected  to  reproduce  the  undulations  of  the  Argand  diagram  of  Q^tv) 
in  Fig  2,  it  would  surely  be  capable  of  representing  the  region  0  <  v  <  1  to  good 


*  This  might  not  be  the  case  for  motion  which  ultimately  diverges,  eg  q(r)  *  t 
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accuracy  and  of  continuing  without  undulations  through  1  <  v  <  3  to  give  usable  results 
for  active  control  studies.  Some  evidence  to  this  effect  is  provided  in  columns  (c)  and 
(d)  of  Table  6,  since  the  former  largely  excludes  the  undulations.  The  present  calcula¬ 
tions  stand  as  a  yard-stick  for  various  approximations  discussed  in  Ref  8. 

5  INVERSE  CALCULATIONS 

In  general,  as  discussed  in  Refs  4,  8  and  14,  the  time-dependent  aerodynamic  forces 
occur  as  terms  in  the  equations  of  an  aeroelastic  system;  for  simplicity,  the  present 
application  is  confined  to  the  aerodynamics  without  any  modelling  of  the  aircraft  response. 
In  practice,  an  active  control  is  subject  to  a  sensor,  such  as  an  accelerometer,  attached 
to  the  aircraft  structure;  within  a  purely  aerodynamic  framework  no  attempt  is  made  to 
simulate  a  sensor  or  any  associated  control  laws.  Instead,  perfect  knowledge  of  the  aero¬ 
dynamic  input  is  assumed,  and  in  section  5.2  this  amounts  to  the  growth  of  lift  due  tc  a 
step  gust. 

In  section  5.1  we  address  the  question  of  determining  the  motion  of  a  control 
surface  to  produce  a  given  time-dependent  aerodynamic  force.  Equation  (2-21) 


dq.(t) 


lj<0>/ 


f  dW 


Fl(T‘T0)dT0 


<V.T> 


is  to  be  solved  for  qj (t)  ,  given  Q^(t)  ,  the  velocity  factor  Fj(o)  from  equation 
(2-22),  and  the  coefficients  Q!j(0)  and  Q'M")  .  The  program,  mentioned  at  the  end  of 
section  3.3,  provides  numerical  values  of  Fj(o)  for  selected  o  up  to  and  including 
some  arbitrary  large  value  Oj  ,  beyond  which  equation  (3-30)  is  used  to  give 

Q'(0)Q!  (0)(o-o  )  |F  (o  )  -  l{(K 

F  (o)  -  1+-^ - ^ - - - L  +  — ! — L_ - L  .  •  (5-2) 

1  8uQ!  j  (0)oJ  <T 


5.1  Alternative  approaches 


One  scheme  of  calculation  is  to  iterate  by  taking  equation  (5-1)  in  the  form 


dq.(T) 

Q!.(0),.W  .  QV.(-)  -L- 


r  dcl-<Tn> 

+  Qij(0)J  ^F1(T“T0)  ■  'K  -  Qi(T) 


and  by  evaluating  qfn+^(t)  from  a  previous  approximation  q^(x)  from  the 
J  J 

relationship 


qjn+0(T) 


<Mt)  QV.(»)  dq^n)(T)  ?dq-n)(T  ) 

Q! .  (0)  ~  QMO)  "^dr  /  JdTfl  iFl(T"V  "  ^dT0 

ij  ij  n  0 


There  is  the  disadvantage  that  each  iteration  involves  the  numerical  differentiation  of 
q^(t)  ,  and  the  procedure  is  often  unstable. 


Instead,  greater  success  is  obtained  from  the  linear  differential  equation 


Q!.(0),'”  °(T>  *  QVj (*)  —  -  Qi(T)  -  q:.(0> ]  - 

A  ” 


1 


1  }dr 


0 

(5-5] 


which  gives 
.  (n+ 1 ) 


qj-  '(t) 


exp[-  Qj {(0)t/QV .(“)] 


ill  ;  _  jj. 
oV  W 


ij 


[  fn  (r*}  -  »YnL 

¥?'’l 

1  Pi' T  '  R  J]  exPI 

o  L 

(5-6] 


where 


,(n) 


J  dqjn)(r  ) 

(x)  -  Q [5(0)J  (f^t-Tq)  -  l}dxc 


(5-7] 


The  iterations  can  be  started  by  taking  qj^  (t)  =  0  and  hence  R^  (t)  =  0  .  Then  by 


equations  (5-6)  and  (5-5)  respectively 


exp  I-  Q 


1.(0)t/QV.(-)]  ?  f*Q! -(0)T*“| 

frw— L JQi(Tl>  •1",-%7=rldT' 

Q  L  1J  -J 


(5-8] 


and 


dqj°(T) 


.0), 


dx 


Q!.(0)q:,;{T)  Q.(x) 

'  Q$/->  +  * 


(5-9] 


Then  equations  (5-7)  to  (5-9)  determine  R^(x)  ,  equation  (5-6)  and  R^  (x)  determine 


.(2), 


qj  (t)  ,  and  its  derivative  is  obtained  by  differentiating  equation  (5-6)  to  give 


dqjn*°(T) 


Q] • (0)q«n+1) (x)  Q.(x)  -  R(n)(x) 


dx 


'Wr 


(5-io: 


with  n»l  .  This  procedure  has  been  used  successfully  for  the  lift  (i *  1)  and  the  modes 
j  *  3  and  5  corresponding  to  the  trailing-edge  and  all-moving-tip  controls,  when  the 
exponent  outside  the  integral  in  equation  (5-6)  is  negative. 

For  the  leading-edge  control  QP(O)  and  QVj(»)  are  of  opposite  sign  and 
numerical  errors,  however  small,  magnify  from  iteration  to  iteration.  This  poses  a 
problem  which  cannot  be  overcome,  for  example,  by  subtracting  2QV ^ («)dqj (x)/dx  from 
both  sides  of  equation  (5-3)  and  allowing  for  the  extra  term  on  the  right  hand  side  in 
the  definition  of  R^  (x)  .  The  difficulty  is  apparently  the  basic  one  that  there 
exists  a  divergent  motion  of  the  leading-edge  control,  which  produces  zero  lift.  There 
are  indications  for  very  high  frequency  from  piston  theory,  considered  in  section  3.2, 
that  the  ultimate  value  of  the  generalized  force  coefficient. 
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q|4(«»)  +  ivQy4(«)  -  -  0.977  +  0.0514  iv  (5-11) 

from  Table  3,  vanishes  for  Che  rapidly  divergent  motion  corresponding  roughly  to 

iv  *  19.0  .  In  this  unusual  situation  any  iterative  method  is  doomed  to  failure  and  a 

different  approach  must  be  adopted. 

We  therefore  take  q4(T)  as  a  series  of  exponential  terms 


q4<T) 


(5-12) 


vaere  the  factors  8^  >  0  are  suitably  chosen  and  the  K  coefficients  are  determined 
so  as  to  satisfy  equation  (5-3)  to  a  sufficient  approximation.  A  good  procedure  is  to 
evaluate  the  left  hand  side  of  equation  (5-3)  for,  say,  2K  values  of  t  in  the  range 
of  interest  and  then  to  solve  the  set  of  linear  simultaneous  equations  that  satisfy  the 
initial  condition  on  q4(0)  and  minimize  the  root  mean  square  error  in  the  quantities 
Q.(t)  . 

5.2  Neutralization  of  gust  force 

We  first  calculate  the  lift  on  the  wing  of  Fig  1  as  it  enters  a  step  gust  at 

M  *  0.8  in  the  direction  normal  to  the  gust  front.  As  discussed  in  section  2.2  and 

reproduced  from  Fig  2b  in  Fig  17,  the  lift  L(o)/L(»)  or  L/Lj  grows  more  slowly  than 
that  due  to  a  step  change  in  the  trailing-edge  control  angle,  which  also  involves 
impulsive  lift  at  o  ■  0  .  Clearly  a  less  violent  change  in  control  angle  is  needed  to 

neutralize  the  growth  of  lift  from  the  step  gust.  The  iterative  calculation,  based  on 

equations  (5-5)  to  (5-10),  gives  the  curve  of  5/5j  with  alternating  long  and  short 

dashes  in  Fig  17.  The  most  striking  feature  of  this  curve  is  that  it  attains  98Z  of  its 

final  value  and  remains  within  ±2.7,  of  it  after  just  four  mean  chords  of  travel,  at  which 
stage  the  lift  from  the  step  gust  is  only  approaching  80Z. 

The  same  method  of  calculation  is  used  for  the  all-moving-tip  control,  but  for  the 
leading-edge  control  satisfactory  results  are  obtained  by  means  of  equation  (5-12)  with 
K  »  4  .  Fig  18  compares  the  neutralizing  curves  of  5/5  j  for  the  three  types  of  control 
surface,  which  show  the  inverse  characteristics  of  their  velocity  factors  Fj(o)  from 
Fig  5.  The  overshoot  in  lift  from  the  leading-edge  control,  typified  by  Fj(o)  >  1  for 

o  <  0.5  ,  leads  to  the  less  rapid  demand  for  5  in  Fig  18,  while  the  all-moving-tip 

shows  intermediate  characteristics,  just  as  it  has  done  in  Fig  5  and  Table  10.  The 
relative  magnitudes  of  5j  for  the  three  controls  in  Fig  18  show  the  high  ratio  10.8 
for  the  leading-edge  control,  which  is  much  less  effective  than  the  trailing-edge  control 
in  producing  lift.  For  each  control  6  remains  within  ±\7,  of  its  final  value  5j  after 
a  gust  penetration  of  only  5.4c,  although  the  growth  of  lift  due  to  the  gust  takes  at 
least  three  times  as  long  to  do  so.  This  remarkable  result  is  peculiar  to  lift  because  of 
the  invariant  asymptotic  behaviour  found  in  equations  (2-17)  and  (3-5).  Nevertheless, 
the  rapid  control  rates  implied  in  Fig  18  are  indicative  of  the  targets  that  can  be  set 
in  active  control  technology.  Such  applications  require  the  full  treatment  of  flow 
history. 
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6  CONCLUDING  REMARKS 
6.1  Recapitulation 

A  double  integral  £or  the  time-dependent  aerodynamic  force  due  to  arbitrary  wing 
or  control-surface  motion  in  equation  (2-19)  forms  the  basis  of  calculations  in  either  the 
time  or  the  frequency  domain.  In  the  time  domain  various  expressions  involving  different 
hereditary  functions  are  derived,  in  particular,  equation  (2-34)  containing  two  quasi¬ 
steady  terms  and  a  'history  term'  involving  the  history  function  H(o)  .  A  program  in 
ICL  1900  Extended  Fortran  has  been  written  to  calculate  the  hereditary  functions  as 
accurately  as  possible.  The  hereditary  functions  for  a  high-aspect-ratio  wing  at  M  ■*  0.8 
with  various  control  surfaces  (Fig  1)  are  used  to  determine  aerodynamic  forces  due  to 
transient  deployment  of  these  controls  in  the  form  of  equation  (4-1).  For  this  particular 
form  of  control  deployment,  analytical  integration  with  respect  to  the  time  variable  leads 
to  the  expression  in  equation  (4-26) ,  which  has  been  evaluated  via  its  Fourier  transform 
in  the  frequency  domain  to  give  good  agreement  with  the  time-dependent  lift  calculated 
from  the  hereditary  functions  in  the  time  domain. 

Special  attention  has  been  paid  to  the  asymptotic  behaviour  of  the  various  functions. 
The  evaluation  of  each  hereditary  function  takes  account  of  the  asymptotic  expansions 
of  the  generalized  forces  Q.j(v)  ^or  small  and  large  frequency  in  equations  (3-1),  (3-2), 
(3-16)  and  (3-17).  The  hereditary  functions  Fn(o)  and  H(o)  are  discussed  in  relation 
to  their  asymptotic  expansions  for  small  and  large  time  in  equations  (3-24)  to  (3-33). 

The  contribution  of  the  history  term  to  lift  due  to  the  trailing-edge  control  is  examined 
in  relation  to  the  limiting  expressions  in  equations  (4-7)  and  (4-10)  for  very  rapid  and 
very  slow  control  rates.  Having  a  well-defined  asymptotic  expansion  for  small  o  ,  the 
history  function  H(o)  is  considered  from  the  standpoint  of  empirical  correction.  The 
factor  <(o)  in  equation  (3-36)  or  (3-37)  is  proposed,  but  there  is  a  simpler  alternative 
in  equation  (3-39). 

Transient  deployment  of  the  controls  according  to  equation  (4-1)  with  variable 
duration  Tj  is  considered,  and  calculations  are  made  of  the  time-dependent  lift, 
pitching  moment  or  hinge  moment  for  pitching  or  control-surface  motion.  The  relative 
importance  of  the  quasi-steady  and  history  terms  is  analysed  in  Tables  7  to  9.  Analysis 
is  also  carried  out  on  the  direct  use  of  oscillatory  aerodynamic  forces  in  approximating 
to  the  time-dependent  lift,  as  in  equation  (4-22).  This  has  considerable  qualitative 
merit,  but  the  quantitative  predictions  can  seriously  exaggerate  the  results. 

The  Fourier  transform  results  for  time-dependent  lift  in  Table  11  confirm  the 
minimal  consequences  of  neglecting  the  logarithmic  term  at  small  frequency  in  equation 
(3-2)  and  show  that,  important  though  it  is  for  the  hereditary  functions,  the  asymptotic 
expansion  for  large  frequency  can  be  ignored  altogether  by  truncating  the  frequency  range 
in  equation  (4-26).  Roughly,  in  the  frequency  domain,  an  upper  limit  of  |v|  «  v  ■  12/Tj 
appears  to  be  satisfactory  (section  4.4). 

Methods  of  inverting  the  equations  to  determine  the  motion  that  will  give  a  required 
aerodynamic  force  are  discussed  in  section  5.1.  The  iterative  scheme  in  equations  (5-5) 
to  (5-10)  appears  to  work  unless  the  coefficients  Q..(v)  imply  a  zero  in  the  lower  half 
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of  the  complex  frequency  plane.  When  this  seems  to  happen  for  the  particular  case  of  lift 
due  to  the  leading-edge  control,  a  more  crude  approach  of  the  collocation  type  has  been 
adopted.  The  illustrative  calculations  of  control-surface  motion  to  neutralize  the 
growth  of  lift  on  entering  a  step  gust  show  different  normalized  control  rates  for  the 
three  types  of  control,  which  should  be  considered  in  relation  to  their  lift  effectiveness 
(Fig  18). 

6.2  General  conclusions 

(l  From  the  consideration  of  both  small  and  large  frequency,  it  is  preferable  to 
calculate  the  hereditary  functions  from  integrals  involving  the  in-quadrature  force 
coefficients  Q'.'.(v)  than  from  alternative  integrals  with  the  in-phase  coefficients 
Q^v)  . 

(2)  When  the  time-dependent  lift  is  expressed  in  terms  of  quasi-steady  displacement  and 
rate  terms  and  the  additional  history  term,  the  history  term  is  likely  to  be  important 
whenever  the  rate  term  is  important. 

(3)  In  calculations  of  time-dependent  hinge  moment  the  two  quasi-  'teady  terms  without 
the  history  term  seem  to  cover  all  practical  needs. 

(4)  During  transient  deployment  of  a  trailing-edge  control  the  streamwise  centre  of 
pressure  of  the  incremental  loading  moves  rapidly  from  an  axis  intersecting  the  control 
to  one  slightly  upstream  of  the  wing  aerodynamic  centre  before  tending  asymptotically  to 
the  aerodynamic  centre  as  the  incremental  lift  disappears. 

(5)  The  time-dependent  lift  due  to  rapid  control-surface  motion  shows  sluggish  behaviour 
in  the  case  of  a  trailing-edge  control,  but  over-reacts  to  a  leading-edge  control.  An 
all-moving  tip  has  intermediate  characteristics  and  produces  lift  that  sometimes  leads  and 
sometimes  lags  the  control  displacement  slightly. 

(6)  The  steady  two-dimensional  chordwise  loadings  for  leading-edge  and  trailing-edge 
controls  provide  a  qualitative  explanation  of  their  time-dependent  lift  characteristics 
on  the  basis  that  the  growth  of  loading  more  remote  from  the  control  surface  is  subject 
to  the  major  delay. 

(7)  For  many  applications  of  active  controls  it  would  be  unnecessary  to  take  into 
account  frequencies  beyond  the  flutter  range.  For  this  approximation  a  much  simpler 
aerodynamic  modelling  would  become  practicable. 

(8)  The  calculated  control-surface  motion  required  to  neutralize  the  growth  of  lift 
on  entering  a  step  gust  stabilizes  within  about  one  third  of  the  corresponding  distance 
for  the  growth  of  lift  itself. 
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FORMULATION  OF  TIME-DEPENDENT  FORCES 


Consider  a  flat-plate  wing  displaced  upwards  from  the  plane  z  ■  0  according  to 
the  equation 

z(x,y,t)  *  Zj (x,y)q(x)  (/ 

where  x  *  Ut/c  is  non-dimensional  time.  The  motion  starts  at  x  *  0  ,  so  that 
q Ct)  *  0  for  x  <  0  and  thereafter 

J  dq(T0) 

q(T)  «  j  “dT-^o  ‘  (* 

n  w 


Remembering  that 


CO 

f  iilLildv 

J  v 


>■  ±  Jn  according  as  x  iO  , 


we  may  write 


!  d'1<,o)  i  i  f  sinp<VT)l  , 

J  —  2-7  J  - = - 4,>dT' 

n  w  a 


,  J  r  dq(x0)  sin[v(TQ-T)] 

"  ¥  J  j  “d^  =  dv 


dvdxQ  . 


The  integrand,  regarded  as  a  function  of  x  ,  corresponds  to  the  sinusoidal  wing  motion 


,  s  Ji  ,  ,  dq(V  i'5(T"T0)l 

Z()(x,y,t)  «  fl^-Zj(x,y)  __e  | 


which  provides  the  generalized  force  coefficient  in  mode  i 


[v  dx0  iQij 


(v)  +  ivQVj (v)}e 


iv(x-xn)l 


—  (=Q!.(v)  si 
i  I  v 


sin  [u(x  -  Xq)]  +  QVj(v)  coslv(x-xQ)] 


Then  the  generalized  force  coefficient  corresponding  to  the  wing  motion  given  by  equations 
(A-l)  and  (A-4)  is 


£RBCtGgjg 


Q*(f)  *  iq(”)Q-j(0)  --  j  J  Q0(x)dvdxQ  . 

~  - —  6  o 


(A-7)  , 


num> 
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It  is  important  to  note  that,  since  the  expression  for  q(x)  in  equation  (A-4)  is 
identically  aero  for  t  <  0  and  covers  the  whole  range  -®  <  t  <  «  ,  so  does  the  quantity 
Q.(x)  defined  by  equations  (A-6)  and  (A-7). 

With  a  little  rearrangement  of  terms  and  the  removal  of  the  subscripts  i  and  ij  , 
we  may  write 

~  dq<Trt) 


f  if  l 

Q(t)  -  J  — — -  |Q*  (0)  +-  J  rQ'(v)  sin(v(x -xQ)]dv 

A  0  rt 


dxr 


oa  am 

.  t  t  dq(ij 

J  J  ~ dr — cos(v(r -TQ)]dvdT0 


0  0 


1  f  f  dq(Tn> 

7Q,,(")J  ]  “d77“  cos[v(T-T0)]dvdT0 

A  A  V 


(A-8) 


0  0 


Given  that  dq(xQ)/dx0  *  0  at  xQ  -  0  and  xQ  -  «  ,  the  factor  of  Q"(»)  in  the  last 
term  of  equation  (A-8)  simplifies  on  integration  by  parts  with  respect  to  xQ  .  Hence 


I  f  f  dq(v 

”  J  J  dTo 

A  A  W 


j  JT  d2q(xQ)  sinlv(x-xQ)] 

CCslv(T-T0)]dvdT0  -  -  J  |  ——ji - = -  dvdT0 


0  0 


0  o  dTo 


1  (  d2q(T0)  .  1  f  d  . 

2  J  2  dT0  "  2  I  ,2  dT 

0  dT0  -  dT" 


00  2 

d  q<*n> 


dq(x) 
dx  * 


if  x  >  0,  (A-9) 


and  vanishes  identically  if  x  <  0  .  Since  Q(x)  must  vanish  identically  when  x  <  0  , 
equation  (A-8)  yields  the  identity 


^  u  n  ' 


Qf(v)  sin  [v(Tq  -  t)]  +  lQu(v)  -Q"(»)]  cosIvCTq-t)] 


dv 


dT0  ’  0 


whatever  the  motion  q(x^)  .  Therefore  for  o  >  0 


00  00 

iQ' (0)  -  —  f  rQ'(v)  sin  vo  dv  +  —  [  lQ"(v)  -  Q"(®)]  cos  vo  dv  «  0  . 

*  J  v  n  J 


(a-io; 


This  identity  allows  the  upper  limits  of  integration  with  respect  to  x^  to  be  changed 
from  ®  to  x  in  the  first  two  terms  of  equation  (A-8),  which  simply  expresses  that 
future  motion  cannot  influence  the  present.  Thus  by  equations  (A-8)  and  (A-9) 
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Q(t)  -  Q"<-) 


J 


CD 

iQ*  (0)  +  j  j  rQ'(v)  sin{v(r -T0))dv 


dt. 


.  J  ?  dq(T.) 

+  n  J  J  — r~  "  Q"(“^  co.|v(t-t0)1 

0  0  0 


dvdtQ  .  (A-!!) 


The  relationship  in  equation  (A-10)  also  allows  the  elinination  of  either  Q'(v)  or 
Q"(v)  frost  the  double  integrals  for  Q(t)  when  r  >  0  ,  so  that  we  obtain  alternative 
expressions 


Q(t)  -  Q"(“)  Q'(0)q(T)  +|  j  J  jQ"(v)  -  Q"(»)f  cos[v(T-T0)]dvdT0 

0  0  0 


(A- 12) 


and 


0M  -  0"<-)  ^  *  T  / /  TT- ^  * i»Bh-T0)ldWT#  .  (*-13) 

0  0  0 
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EVALUATION  OF  HEREDITARY  FUNCTIONS 

Three  Integrals 

The  hereditary  functions  in  equations  (3-23)  involve  the  following  integrals 


where 


v°) 


I,(o) 


I2(o) 


J f(v)v  sin  vo  dv 
0 

ao 

J f(v)  cos  vo  dv 
0 

00 

[  f(v)  lisis.  dv 
J  V 

0 


> 


J 


f(v) 


QV.(v)  -  QV.<») 

qttto r— 


(B-l) 


(B-2) 


Each  integral  is  evaluated  in  three  parts  by  subdividing  the  range  of  integration  into 
0<v<v,  v<v<v  and  v  <  v  <  «  .  Different  integration  procedures  are  used  in 

XX  u  u 

the  three  regions. 


0  <  v  <  v„ 


As  explained  in  section  2.1,  QVj(v)  is  given  by  equation  (3-2)  so  that  equation 


(B-2)  becomes 


f(v) 


2  _2 

B0  “  Q'i j  +  B'j'v  +  B2v  log  v  +  B'jV 

- U - q]77ol - 


(B-3) 


’U 


Typically  we  take  v^  *  0.08,  subdivide  the  range  into  16  equal  parts  and  evaluate  the 
contributions  to  the  integrals  (B-l)  by  Simpson's  Rule,  for  example 


0.08 
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f(v)  cos  vo  dv 
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j  f(0)  + 
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3  V  200  J 


cos 


(2k-  l)o 
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V.  <  V  <  V 
&  u 

All  the  available  values  of  QV.(v)  from  lifting-surface  theory,  excluding  v  *  0  , 
are  used  to  fit  a  standard  cubic  spline.  Typically  the  region  0.08  <  v  <  6.00  is  sub¬ 
divided  into  2n  -  148  equal  parts  and  the  integrands  of  (B-l)  are  evaluated  without  the 
trigonometrical  factor  sin  vo  or  cos  vo  .  They  are  then  represented  by  quadratics  over 
successive  double  intervals 


w4  +  (2k-2)h  <  v  <  v  +  2kh 


(h  »  »  k  *  1 ,2 . n^  (B-5) 


as  for  Simpson’s  Rule,  but  in  the  subsequent  integration  the  trigonometrical  factor  is 
treated  analytically.  The  procedure  is  known  as  Filon's  Rule  and  is  used  in  preference 
to  Simpson's  Rule  to  ensure  the  necessary  accuracy  when  o  is  large  and  the  trigono¬ 
metrical  factor  is  of  small  wavelength.  The  resulting  contributions  to  Ij(o)  and  I^o) 
are  as  follows: 
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where  with  0  *  oh 


2  2 
0+8  sin  0  cos  8-2  sm  8 


2[8(1 +  cos  o)  -  2  sin  8  cos  8] 

.3 


4[sin  8-8  cos  8] 
«3 


Iq(o)  is  treated  in  the  same  way  as  12(0)  with  vf(v)  replacing  f(v)/v  throughout 
equation  (B-7). 

v  <  v  <  • 
u _ 

As  explained  in  section  3.2,  Q'.'j(v)  is  given  by  equation  (3-17)  so  that  equation 

(B-2)  becomes 

f(v)  -  |“[  +  4j/Qij(0>  *  (B"9 


The  integrals  (B-l)  are  evaluated  analytically  in  terms  of  the  function 


00 

f  ^dp  . 

Z  P 


(B-10) 


Successive  integration  yields  the  results 
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q:.(o)  [ f(v)  dv 

ij  J  V 


(  A"a^\  A" 

-  \A"o2 - S3(vuo)  +  — ^3  sin  vuo  +  cos  \>uoj  .  (B-l 3) 


The  quantity  S3(x)  from  equation  (B-10)  is  obtained  from  expansions  in  series  of 
Chebyshev  polynomials  by  the  procedure  described  in  Ref  15.  The  formulae,  due  to  Davies 
and  hitherto  unpublished,  are 
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for  0  <  x  <  7 


and  for  7  <  x  <  » 


S3(x) 


sin  x 

T~ 


x 


r-0 


COS  X 

5"“ 


X 


(B-14) 


(B-15) 


where  the  dash  '  beside  the  summation  sign  indicates  that  the  term  r  ■  0  is  halved. 
The  Chebyshev  polynomials  T^  (z)  are  defined  in  section  3  of  Ref  15  and  are  evaluated 
by  a  recurrence  relationship.  The  coefficients  Ar,  Br  and  Dy  are  given  in  Table  12 
and  ensure  high  accuracy. 


i  or  j 


Table  1 

DEFINITION  OF  FORCE  AND  UPWASH  MODES 


Force  mode 

Upwash  mode 

Force 

z^/c 

Motion 

v./u 

Lift 

1 

Heaving 

iv 

Pitching  moment 

(x-x  )/c 
a 

Pitching 

1  +  iv(x-x  )/c 
a 

Hinge  moment 

(x-xh3)/c 

Trailing-edge  control 

!  +  iv(x-x^3)/c 

Hinge  moment 

X 

1 

m* 

.O 

V-/ 

01 

Leading-edge  control 

1  +  iv(x-xh4)/c 

Hinge  moment 

(x“*h5)/S 

All-moving  tip 

1  +  iv(x-x^3)/c 

Bending  moment 

|y|/« 

Sinusoidal  gust 

exp(ivx/c) 

The  axis  positions  xa»  x^-jCy),  *^4^  atK*  *h5  an<*  t^e  sPanw*se  extent  °*  tae 
control  surfaces  y^  <  y  <  s  are  defined  in  Fig  1.  In  nodes  3,  4  and  5  the  fomulae 
for  z./c  and  v./U  apply  over  the  extent  of  the  control  surface,  but  both  quantities 
are  nominally  zero  over  the  rest  of  the  planform  (y  <  y^  or  x  <  ®°de  3,  or 

x  >  xh/,(y)  in  mode  4). 


Table  2 


FORCE  COEFFICIENTS  FOR  OSCILLATING  T RAILING-EDGE  CONTROL  (M  -  0.8) 


V 

Lift 

Pitching  moment 

Hinge  moment 

-*13 

-0" 

M13 

-«23 

MM 

~Q33 

-0" 

y33 

0 

-2.6491 

-1.1691 

0.01827 

0.00946 

0.05 

1.7625 

-2.4487 

-1.0902 

0.01827 

0.00952 

0.15 

1.6350 

-1.9144 

-0.8954 

0.01826 

0.00964 

0.40 

1.3271 

-0.8774 

1.7035 

-0.5526 

0.01834 

0.00987 

0.60 

1.2257 

-0.4838 

1.6210 

-0.4072 

0.01856 

0.00995 

0.75 

1.2083 

-0.3544 

1.5701 

-0.3232 

0.01881 

0.00993 

1,00 

1.1707 

-0.2967 

1.5220 

-0.2324 

0.01929 

0.00981 

1.20 

1.0927 

-0.2464 

IB 

-0.1987 

0.01974 

0.00979 

1.50 

1.0652 

-0.1638 

BS9 

-0.1555 

0.02082 

C. 00966 

1.80 

1 .0040 

-0.1582 

flKH 

-0.1340 

0.02181 

0.00936 

2.10 

0.9495 

-0.1026 

1.3341 

-0.1043 

0.02316 

0.00913 

2.40 

0.9259 

-0.0977 

1.2891 

-0.0889 

0.02425 

0.00868 

2.70 

0.8691 

-0.0672 

1.2427 

-0.0682 

0.02533 

0.00840 

3.00 

0.8485 

-0.0557 

1.1964 

-0.0545 

0.02616 

0.00800 

3.30 

0.8142 

-0.0382 

1.1580 

-0.0372 

0.02673 

0.00777 

3.60 

0.7931 

-0.0262 

1.1233 

-0.0235 

0.02721 

0.00756 

3.90 

0.7784 

-0.0143 

-0.0102 

0.02767 

0.00745 

4.20 

0.7683 

-0.0049 

1.0819 

+0.0015 

0.02830 

0.00739 

4.50 

0.7673 

•*•0.0023 

1.0738 

0.0104 

0.02922 

0.00732 

4.80 

0.7662 

0.0082 

1 .0700 

0.0182 

0.03024 

0.0C722 

5.10 

0.7711 

0.0117 

1.0706 

0.0229 

0.03135 

0.00709 

5.40 

0.7694 

0.0145 

1.0690 

0.0270 

0.03240 

0.00693 

5.70 

0.7725 

0.0168 

1.0699 

0.0299 

0.03331 

0.00675 

6.00 

0.7735 

0.0177 

1.0712 

0.0317 

0.03406 

0.00658 

00 

0.582 

0.0514 

0.874 

0.0735 

0.0354 

0.00516 

-Bi  - 

1 .8046 

- 

0.7462 

- 

0.00058 

- 

j] 


Table  3 

FORCE  COEFFICIENTS  FOR  OSCILLATING  LEADING-EDGE  CONTROL  (M  -  0.8 


Lift 

Pitching  moment 

Hinge  moment 

V 

-0" 

-«24 

-0" 

^24 

-o1 

-0" 

W44 

0 

DBi 

-0.0776 

0.0220 

0.0847 

-0.05667 

0.00312 

0.05 

HU 

-0.0609 

0.0214 

0.0915 

-0.05666 

0.00294 

0.15 

0.1558 

-0.0228 

0.0190 

0.1065 

-0.05661 

0.00252 

wmm 

0.1357 

+0.0406 

0.0174 

0.1284 

-0.05672 

0.00181 

IKS 

0.1243 

0.0714 

0.0212 

0.1318 

-0.05720 

0.00166 

IS 

0.1209 

0.0923 

0.0236 

0.1320 

-0.05775 

0.00184 

1.00 

0.1372 

0.1130 

0.0284 

0.1379 

-0.05861 

0.002*44 

1.20 

0.1559 

0.1072 

0.0408 

0.1424 

-0.05922 

0.00270 

1.50 

0.1564 

0.1071 

0.0630 

0.1424 

-0.06076 

0.00330 

1.80 

0.1986 

0.1172 

0.0948 

0.1435 

-0.06134 

0.00418 

2.10 

0.2252 

0.  >02? 

0.1291 

0.1681 

0.1381 

-0.06209 

0.00453 

2.40 

0.2505 

0.1080 

0.1365 

-0.06222 

0.00521 

2.70 

0.2988 

0.1009 

0.2126 

0.1309 

-0.06196 

0.00533 

3.00 

0.3333 

0.0969 

0.2604 

0.1268 

-0.06203 

0.00555 

3.30 

0.3782 

0.0921 

0.3137 

0.1206 

-0.06185 

0.00567 

3.60 

0.4244 

0.0867 

0.3688 

0.1140 

-0.06178 

0.00578 

3.90 

0.4711 

0.0807 

0.4252 

0.1070 

-0.06172 

0.00590 

4.20 

0.5206 

0.0746 

0.4847 

0.0994 

-0.06141 

0.00604 

4.50 

0.5700 

0.0678 

0.5444 

0.0914 

-0.06095 

0.00616 

4.80 

0.6197 

0.0607 

0.6046 

0.0825 

-0.06019 

0.00624 

5.10 

0.6665 

0.0533 

0.6609 

0.0734 

-0.05941 

0.00628 

5.40 

0.7112 

0.0460 

0.7162 

0.0643 

-0.05845 

0.00631 

5.70 

0.7544 

0.0386 

0.7678 

0.0549 

-0.05734 

0.00631 

6.00 

0.7928 

0.0309 

0.8146 

0.0454 

-0.05608 

0.00626 

GO 

0.977' 

-0.0514 

0.982 

-0.0477 

-0.0531 

0.00516 

-B 


0.1674 


0.0692 


-0.00165 
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Table  4 

FORCE  COEFFICIENTS  FOR  OSCILLATING  ALL-MOVING  TIP  (M  «  0.8) 


V 

Lift 

Pitching  moment 

Hinge  moment 

-^5 

-Q'i's 

-Q25 

-0" 

W25 

~^55 

-Q55 

0 

1.0368 

-0.9826 

0.9567 

-0.0945 

0.00029 

0.06800 

0.05 

1.0239 

-0.8739 

0.9520 

-0.0514 

0.00027 

0.06800 

0.15 

0.9616 

-0.6036 

0.9306 

+0.0478 

0.00016 

0.06797 

0.40 

0.8128 

-0.0995 

0.8816 

0.2092 

-0.00078 

0.06821 

0.60 

0.7622 

+0.1018 

0.8548 

0.2675 

-0.00208 

0.06890 

0.75 

0.7604 

0.1778 

0.8382 

0.3008 

-0.00322 

0.06946 

1.00 

0.7777 

0.2039 

0.8314 

0.3429 

-0.00570 

0.07029 

1.20 

0.7546 

0.2071 

0.8369 

0.3544 

-0.00834 

0.07152 

1.50 

0.7520 

0.2552 

0.8423 

0.3672 

-0.01123 

0.07378 

1.80 

0.7849 

0.2361 

0.8604 

0.3661 

-0.01485 

0.07550 

2.10 

0.7418 

0.2589 

0.8605 

0.3720 

-0.01587 

0.07841 

2.40 

0.8075 

0.2542 

0.8880 

0.3695 

-0.01660 

0.07940 

2.70 

0.7668 

0.2572 

0.8947 

0.3722 

-0.01471 

0.08196 

3.00 

0.8182 

0.2587 

0.9233 

0.3659 

-0.01170 

0.08193 

3.30 

0.8055 

0.2573 

0.9290 

0.3674 

-0.00788 

0.08357 

3.60 

0.8363 

0.2562 

0.9579 

0.3619 

-0.00144 

0.08285 

3.90 

0.8356 

0.2549 

0.9649 

0.3605 

+0.00241 

0.08306 

4.20 

0.8521 

0.2534 

0.9820 

0.3572 

0.00768 

0.08244 

4.50 

0.8600 

0.2536 

0.9984 

0.3566 

0.01262 

0.08229 

4.80 

0.8764 

0.2512 

1.0165 

0.3528 

0.01713 

0.08148 

5.10 

0.8854 

0.2508 

1.0324 

0.3517 

0.02113 

0.08134 

5.40 

0.9036 

0.2486 

1.0539 

0.3487 

0.02697 

0.08079 

5.70 

0.9140 

0.2465 

1.0720 

0.3458 

0.03208 

0.08023 

6.00 

0.9258 

0.2442 

1.0859  ' 

0.3427 

0.03680 

0.07952 

00 

0.855 

0.2331 

C.986 

0.3271 

0.0111 

0.07557 

-b;  « 

1.0465 

- 

0.4327 

- 

0.00026 

- 
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Table  8 

USEFULNESS  OF  QUASI-STEADY  APPROXIMATIONS  TO  TRANSIENT  LIFT  AND  HINGE  MOMENT 
(Errors  are  given  as  percentages  of  the  peak  steady  value) 


Motion 

Force 

Tl 

One-term 

approximation 

Two-term 

approximation 

Maximum 

error 

Rms 

error 

Maximum 

error 

Rms 

error 

Pitching 

Lift 

5 

35.9 

20.8 

67.1 

41.1 

10 

18.3 

10.1 

32.4 

19.2 

20 

8.4 

4.8 

13.8 

8.0 

40 

3.7 

2.3 

4.9 

2.9 

80 

1.8 

1.2 

1.5 

0.9 

Trai ling-edge 

Lift 

5 

31.2 

18.0 

103.8 

62.9 

control 

10 

23.8 

14.6 

48.5 

28.7 

20 

16.8 

10.9 

19.7 

11.6 

40 

10.4 

7.0 

6.7 

4.0 

Trailing-edge 

Hinge 

5 

37.0 

24.7 

5.4 

3.6 

control 

moment 

10 

18.6 

12.5 

1.5 

0.9 

20 

9.1 

6.2 

0.4 

0.2 

40 

4.5 

3.0 

0.1 

0.1 

Leading-edge 

Lift 

20 

9.3 

5.4 

14.2 

8.4 

control 

40 

3.9 

2.4 

5.0 

3.0 

Leading-edge 

Hinge 

20 

0.7 

0.4 

0.4 

0.2 

control 

moment 

40 

0.4 

0.3 

0.1 

0.1 

All-moving  tip 

Lift 

20 

11.1 

7.1 

16.9 

10.0 

6.5 

4.3 

5.8 

3.5 

Table  9 

ANALYSIS  OF  ROOT  MEAN  SQUARE  ERRORS  WITH  QUASI-STEADY  APPROXIMATIONS 


Motion 

Force 

Equal  with  one 
and  two  terms 

Errors  below  5 7, 

Rms  error 

T1 

With 

ane  term 

With 

two  terms 

Pitching 

Lift 

1.67. 

55 

Tl 

> 

19 

TI 

>  28 

Trailing-edge 

control 

Lift 

107. 

23 

T! 

> 

65 

Tl 

>  36 

Leading-edge 

control 

Lift 

2.0% 

50 

T1 

> 

21 

T1 

>  28 

All-moving  tip 

Lift 

5.27. 

31 

T! 

> 

33 

T1 

>  32 

Trailing-edge 

control 

Hinge 

moment 

- 

- 

T1 

> 

24 

TI 

>  4 

Leading-edge 

control 

Hinge 

moment 

- 

T1 

> 

4 

Tl 

>  2 

Tabic  10 


TRANSIENT  LIFT  DUE  TO  WING  AND  CONTROL-SURFACE  MOTIONS 


0069 
0467 
0.2621 
0.5927 
0.8847 
1 .0000 
0.8847 
0.5927 
0.2621 
0.0467 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 
0.0000 


— 

Values 

of  L/Lq 

Pitching 

Trail ing-edge 
control 

Leading-edge 

control 

All-moving 

tip 

0.0439 

0.0050 

0.0089 

0.0232 

0.1567 

0.0280 

0.1069 

0.0943 

0.4875 

0.1529 

0.5213 

0.3420 

0.7996 

0.3624 

0.9145 

0.6170 

0.9397 

0.5782 

1 .0558 

0.7891 

0.8463 

0.7082 

0.8865 

0.7835 

0.5660 

0.6981 

0.4962 

0.6085 

0.2334 

0.5580 

0.0734 

0.3502 

0.0102 

0.3590 

-0.1735 

0.1362 

-0.0118 

0.1974 

-0.1282 

0.0628 

+0.0532 

0.1251 

+0.0893 

0.0879 

0.0424 

0.0821 

0.1033 

0.0726  ! 

0.0351 

0.0539 

0.0648 

0.0494 

0.0115 

0.0116 

0.0131 

0.0042 

0.0021 

0.0027 

0.0015 

0.0026 

0.0004 

0.0004 

0.0004 

0.0004 

Values  of 


0.0005 

0.0069 

0.0467 

0.2621 

0.5927 

0.8847 

1.0000 

0.8847 

0.5927 

0.2621 

0.0467 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 


Pitching 

Trail ing-edge 
control 

Leading-edge 

control 

All-moving 

tip 

0.0012 

0.0003 

0.0010 

0.0008 

0.0102 

0.0043 

0.0107 

0.0076 

0.0538 

0.0320 

0.0554 

0.0441 

0.2638 

0.2019 

0.2641 

0.2347 

0.5756 

0.4940 

0.5744 

0.5363 

0.8513 

0.7870 

0.8495 

0.8197 

0.9638 

0.9466 

0.9621 

0.9547 

0.8614 

0.8994 

0.8608 

0.8791 

0.5925 

0.6677 

0.5931 

0.6286 

0.2848 

0.3614 

0.2863 

0.3221 

0.0800 

0.1238 

0.0816 

0.1018 

0.0231 

0.0325 

0.0232 

0.0278 

0.0062 

0.0072 

0.0062 

0.0065 

0.0025 

0.0028 

0.0025 

0.0026 

0.0007 

0.0008 

0.0007 

0.0007 

0.0002 

0.0003 

0.0002 

0.0002 

Table  11 


LIFT  DUE  TO  RAPID  DEPLOYMENT  OF  TRAILING-EDGE  CONTROL  (t j  -5) 


Full* 

calculation 

Truncation*  in  frequency  parameter  at: 

r 

v  ■  4.91 
u 

v  *  2.45 
u 

v  *  1.84 
u 

V  *  1.23 
u 

0 

0.0000 

-0.0010 

-0.0066 

-0.0280 

0.0326 

0.1 

0.0006 

+0.0002 

-0.0039 

-0.0203 

0.0517 

0.5 

0.0280 

0.0288 

+0.0325 

+0.0410  j 

0.1411 

1.0 

0.1529 

0.1528 

0.1559 

0.1871 

0.2720 

1.5 

0.3624 

0.3631 

0.3601 

0.3829 

0.4035 

2.0 

0.5782 

0.5787 

0.5777 

0.5686 

0.5113 

2.5 

0.7082 

0.7088 

0.7114 

0.6792 

0.5748 

3.0 

0.6981 

■zia 

0.6983 

0.6755 

0.5821 

3.5 

0.5580 

HSI 

0.5548 

0.5640 

0.5330 

4.0 

0.3590 

B9 

0.3612 

0.3924 

0.4391 

4.5 

0.1974 

0.1983 

0.2033 

0.2248 

0.3205 

5.0 

0.1251 

0.1239 

0.1193 

0.1106 

0.2008 

5.5 

0.0964 

0.0971 

0.0927 

0.0643 

0.1012 

6.0 

0.0821 

0.0815 

0.0859 

0.0666 

0.0354 

7.0 

0.0650 

0.0650 

0.0631 

0.0872 

0.0115 

8.0 

0.0538 

0.0540 

0.0538 

0.0488 

0.0632 

10.0 

0.0396 

0.0386 

0.0366 

0.0497 

0.0500 

12.0 

0.0116 

0.0120 

0.0104 

0.0016 

-0.0087 

16.0 

0.0045 

0.0051 

0.0065 

0.0030 

-0.0098 

20.0 

0.0012 

0.0019 

0.0010 

0.0069 

+0.0064 

Rms  error 

0.5  <  t  <  5.0 

0.17. 

0.4% 

2.3% 

9.6% 

*  Total  force  as  in  Table  7(a). 

*  Fourier  transform  calculations  with  finite  range  of  integration  0  <  v  «  v 


Table  12 


COEFFICIENTS  IN  EQUATIONS  (B-14)  AND 


r 

A 

r 

B 

r 

0 

5.18081405 

-1.81985802 

91 

-0.36052882 

0.08176881 

MM 

0.04113346 

-0.00718178 

B 

-0.00635944 

0.00092653 

1M 

0.00121010 

-0.00015385 

5 

-0.00026811 

0.00003058 

6 

0.00006679 

-0.00000697 

7 

-0.00001827 

0.00000177 

8 

0.00000540 

-0.00000049 

9 

-0.00000170 

0.00000015 

10 

0.00000057 

-0.00000004 

11 

-0.00000020 

0.00000002 

12 

0.00000007 

-0.00000001 

13 

-0.00000003 

- 

14 

0.00000001 

- 

7.1517184311 

2.2011779268 

-0.3109142653 

0.0555982922 

-0.0074123262 

0.0007046314 

-0.0000490713 

0.0000025918 

-0.0000001070 

0.0000000035 

-0.0000000001 
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LIST  OF  SYMBOLS 


A  ,B  ,D 
r’  r’  r 


*K*  R*  T 


coefficients  (p  «  0,1,2)  in  equation  (3-13) 

coefficients  (p  =  0,1,2)  in  equation  (3-16) 

coefficients  (p  *  0,1,2)  in  equation  (3-17) 

coefficients  (r  *  0,1,2,....)  in  Table  12 

coefficients  (p  =  0,1, 2, 3)  in  equation  (3-13) 

coefficients  (p  =  0, 1,2,3)  in  equation  (3-1) 

coefficients  (p  *  0,1, 2, 3)  in  equation  (3-2) 

local  chord  of  wing 

mean  chord  of  wing  (see  Fig  1) 

crank  chord,  root  chord,  tip  chord  (see  Fig  1) 

subscript  denoting  empirical  value 
function  in  equation  (B— 2) 

hereditary  position  factor  in  equations  (2-25) 

hereditary  velocity  factor  in  equation  (2-22)  or  (2-23) 

hereditary  acceleration  factor  in  equation  (2-31) 

planform  rounding  function  in  equation  (2-7) 
subdivided  interval  in  equation  (B-5) 
history  function  in  equation  (2-35) 

/-T;  integer  or  subscript  denoting  force  mode  (see  Table  1) 
integral  in  equation  (4-24)  or  (4-25) 
integrals  (p  =  0,1,2)  in  equations  (B-l) 

integer  or  subscript  denoting  upwash  mode  (see  Table  1) 
hereditary  function  in  equation  (4-12) 
variable  integer 

Kiissner  kernel  function  in  equation  (2-3) 
number  of  terms  of  series  in  equation  (5-12) 
complex  wing  loading  in  equation  (2-4) 

instantaneous  lift  force  in  equations  (2-13)  and  (4-22) 
peak  value  of  quasi-steady  lift  force 

final  steady  lift  force 

number  of  spanwise  terms 
Mach  number  of  stream 
instantaneous  pitching  moment 
superscript  denoting  nth  iteration 
number  of  chordwise  terms 
integers  (i  *  1,2,...,N)  (see  Ref  5) 
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LIST  OF  SYMBOLS  (continued) 

time-dependent  factor  for  mode  j  in  equations  (2-18),  (4-1)  and  (A-l) 

elementary  generalized  force  coefficient  in  equation  (A-6) 

time-dependent  force  coefficient  in  mode  i 

complex  generalized  force  coefficient  in  equation  (2-6) 

real  part  of 

imaginary  part  of  Q.,j/v 
quantities  in  equation  (4-3) 

iterative  integral  in  equation  (5-7) 

real  part  of 
semi-span  of  wing 
integral  in  equation  (B-10) 

time 

duration  of  wing  or  control-surface  motion 
Chebyshev  polynomials  (r  *  0,1,2,,...) 
speed  of  stream 

upwash  angle  in  mode  j  (see  equation  (2-2)  and  Table  1) 

streamwise  coordinate  (see  Fig  1) 
location  of  pitching  axis  (see  Fig  1) 

local  x-coordinate  of  hinge  line  (j  ■  3,4,5  in  Fig  1) 

x-coordinate  of  local  leading  edge 

location  of  tip  leading  edge  (see  Fig  1) 

x-coordinate  of  local  trailing  edge 

spanwise  coordinate  (see  Fig  1) 

y-coordinate  of  inboard  end  of  control  surface  (i  *  3,4,5) 

y-coordinate  of  crank  station  (see  Fig  1) 

upward  displacement  from  planform 

elementary  sinusoidal  wing  motion  in  equation  (A-5) 

force  mode  in  Table  1,  displacement  in  equation  (2-1) 

functions  in  equations  (B-8) 
coefficients  and  factors  in  equation  (5-12) 

instantaneous  control  deflection  angle  in  streamwise  plane  (radians) 
peak  value  of  6  (radians  unless  otherwise  stated) 

final  control  deflection  angle 
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LIST  OF  SYMBOLS  (concluded) 

2 

60^  lift  per  unit  area/JpU 

AL  contribution  to  L  from  integral  in  equation  (2-34) 

6  instantaneous  angle  of  pitch  with  heave  constraint 

6q  peak  value  of  6 

k,k  empirical  correction  factors 

X  planform  rounding  variable  in  equations  (2-8)  and  (2-9) 

v  frequency  parameter,  wc/U 

v  ,v  lower  and  upper  values  subdividing  the  range  0  <  v  <  ® 

u 

p  density  of  stream 

o  distance  of  travel  in  mean  chords,  Ut/c 

Oj  value  of  a  above  which  equation  (5-2)  is  used  for  Fj 

om  value  of  o  ,  if  any,  where  H(o)  is  a  maximum 

t  non-dimensional  time,  Ut/c 

t  j  non-dimensional  duration  of  motion,  Utj/c 

T2  value  of  x  when  L/Lq  has  fallen  to  0.01 

angular  chordwise  parameters  in  equations  (4-21) 

w  circular  frequency 
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Fig  2a&b  Normalized  lift  from  oscillatory  and  step  motions 


Fig  4a  Three  hereditary  function! 


Fig  5  Velocity  factors  for  lift  in  various  motions 
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Fig  8  Time-dependant  forces  for  trailing-edge  centre!  motion 


Fig  9  Contributions  to  lift  from  trailing -adga  control 
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